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Preface 1 



In my tenth year at the Institute, I dedicate this book to the BCIT community 



The primary purpose of writing a book and distributing it free-of-charge is to extend my gratitude 
to BCIT 2 . I am particularly thrilled to do it with this textbook because it is a product of many 
learning opportunities BCIT has offered me over a period of several years. What follows is a brief 
background on how this book came to be. 

My post-secondary teaching career began on 22 January 2001 at the Pacific Marine Training Cam- 
pus of BCIT when I logged on to a Unix workstation to instruct in the Propulsion Plant Simulator. 
That has been a major milestone in many ways in my professional life. While learning inner 
workings of Unix operating system (OS), I also made a discovery and that discovery profoundly 
changed my view on how I thought the world operated. The discovery was the GNU/Linux OS 
and open source software (OSS) movement through several books, most notably Just for Fun: The 
Story of an Accidental Revolutionary 3 and The Cathedral and the Bazaar 4 . I was convinced that 
the collective power of connected individuals around the world and the global infrastructure of the 
Internet had the potential to change the ways the world functioned. 

In the last 10 years, BCIT has allowed me to study various subjects through its Professional De- 
velopment (PD) programs for which I am very grateful. I learned a great deal in PD courses and 
in one of the recent ones, I had two deja vu moments similar to my discovery of OSS movement. 
The first one occurred when I began reading The Wealth of Networks 5 and the second one when 
I found about Connexions 6 . The former was a confirmation of my 10-year old discovery and the 
latter is what I am using to write this book. Connexions is a web-based curricular content author- 
ing and publishing technology that I believe has a growing potential for writing and distributing 
free-of-charge learning materials. 

Thus, motivation for this book stems from the notions that were generated by the OSS movement. 



'This content is available online at <http://cnx.org/content/m4 1458/1. 6/>. 
2 http://www.bcit.ca/ 

3 Just for Fun: The Story of an Accidental Revolutionary by L. Torvalds and D. Diamond, New York: HarperCollins 
Publishers. ©2001 

4 The Cathedral and the Bazaar by E. S. Raymond, Sebastopol: O'Reilly Media. ©1999 
5 The Wealth of Networks by Y. Benkler, New Haven: Yale University Press. ©2006 



6 http://cnx.org/ 



The book was written to pay a small token of appreciation to BCIT and I hope it will be a contri- 
bution to the open educational resources repository. 

Serhat Beyenir 

North Vancouver, B.C. 

25 October 2011 



Study Guide 



MATLAB, a sub-course of Computer Technology 1 and this text are specifically designed for 
students with no programming experience. However, students are expected to be proficient in First 
Year Mathematics and Sciences and access to good reference books are highly recommended. I 
also assume that students have a working knowledge of the Mac OS X or Microsoft Windows 
operating systems. 

The strategic goal of the course and book is to provide learners with an appreciation for the role 
computation plays in solving engineering problems. The MATLAB specific skills that I would like 
students to acquire are as follows: 

• Write scripts to solve engineering problems including interpolation, numerical integration 
and regression analysis, 

• Plot graphs to visualize, analyze and present numerical data, 

• Publish reports. 

The best way to learn about engineering computation is to actually do it. We will therefore solve 
many engineering problems mainly using a recent version of MATLAB in this book. Since the 
primary focus is engineering computation, we will concentrate on the mathematical solutions and, 
to a limited extent, the graphical user interface (GUI) features of MATLAB. 

Learning a new skill, especially a computer program, can be an overwhelming experience. To 
make the best of this process, students are encouraged to observe the following guidelines that 
have proven to work well: 

• Plan to study 2 hours outside of class for every hour inside of class, 

• Practice, practice, practice: As the old saying goes, practice makes one perfect or perhaps 
we should modify that statement: Good practice makes one perfect, 

• Buddy system: Study with a classmate. Helping one another drastically improves your 
understanding of the material. Particularly, students are advised to work the problem sets in 
this fashion, 

• Muddy points: Make a note of muddy points as they may occur during lectures and email 
your notes to me. I will address those issues at the beginning of the next class, 

• Open book exam: Do not try to memorize commands, functions or their syntax but learn 
where and how to find that information. Through many exercises and problem sets you will 



7 This content is available online at <http://cnx.org/content/m4 1459/1 .2/>. 



have solved by the end of the course, most computational routines will become second nature 
to you. The exam is open book, so keep your learning materials and m-files well organized. 



Chapter 1 
Introduction 



1.1 What is MATLAB? 



V V WTUSR 




MATLAB stands for MATrix LABoratory (see wikipedia 2 ) and is a commercial software appli- 
cation written by The MathWorks, Inc. 3 When you first use MATLAB, you can think of it as 



lr rhis content is available online at <http://cnx.org/content/m4 1403/1. 2/>. 
2 http://en.wikipedia.org/wiki/MATLAB 

3v 



http://www.mathworks.com/ 



6 CHAPTER 1. INTRODUCTION 

a glorified calculator allowing you to perform engineering calculations and plot data. However, 
MATLAB is more than an advanced scientific calculator, for example MATLAB's sophisticated 
numerical computation environment also allows us to analyze data, simulate engineering systems, 
document and share our code with others. 

1.1.1 Why Use MATLAB? 

MATLAB has become a defacto standard in many fields of engineering and science. Even a ca- 
sual exploration of MATLAB should unveil its computational power however a closer look at 
MATLAB's graphics and data analysis tools as well as interaction with other applications and 
programing languages prove why MATLAB is a very strong application for technical computing. 

The standard MATLAB installation includes graphics features to visualize engineering and scien- 
tific data in 2-D and 3-D plots. We can interactivity build graphs and generate MATLAB command 
output that can be saved for use in the future. The saved-instructions can be called again with dif- 
ferent data set to build new plots. The plots created with MATLAB can be exported in various file 
formats (e.g. .jpg, .png) to embed in Microsoft Word documents or PowerPoint slideshows. 

MATLAB also contains interactive tools to explore and analyze data. For example, we can visu- 
alize data with one of the many plotting routines, zoom in to plots to take measurements, perform 
statistical calculations, fit curves to data and evaluate the obtained expression for a desired value. 

MATLAB interacts with other applications (e.g. Microsoft Excel) and can be called from C code, 
C++ or Fortran programming language. 

1.1.2 Running MATLAB 

To use MATLAB, it must be installed on your computer and you can start it just like you start any 
application on your system or you must have access to a network where it is available. 

In POWR 3307, we will use MATLAB by accessing the BCIT network. The network access 
is platform independent, that is, we can run MATLAB under Mac OS X or Microsoft Windows 
operating systems through a web browser. The following links provide instructions on how to 
access and use BCIT's AppsAnywhere service: 

How to access AppsAnywhere with Safari on a Macintosh Computer 4 

How to open and save files in AppsAnywhere when logging in from a Macintosh 5 

How to access AppsAnywhere using Firefox 6 

How to open and save files in AppsAnywhere when logging in from Windows 7 



4 https://helpdesk.bcit.ca/fsr/sr/appsany where/750. html 
5 https://helpdesk.bcit.ca/fsr/sr/appsany where/807. html 
6 https://helpdesk.bcit.ca/fsr/sr/appsany where/701. html 
7 https://helpdesk.bcit.ca/fsr/sr/appsany where/806. html 



A trial version of MATLAB can be obtained from the mathworks website. 8 

1.1.3 The MATLAB Desktop 

When you start the MATLAB program, it displays the MATLAB desktop. The desktop is a set of 
tools (graphical user interfaces or GUIs) for managing files, variables, and applications associated 
with MATLAB. The first time you start MATLAB, the desktop appears with the default layout, as 
shown in the following illustration. 
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Figure 1.1: The MATLAB Desktop. 



1.1.3.1 Command Window 



The Command Window is where we execute MATLAB commands. We enter statements at the 
Command Window prompt. The prompt can be any one of the following: 



s http://www.mathworks.com/products/matlab/tryit.html 
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Trials indicates that the Command Window is in normal mode and the MATLAB license 
will expire after the trial period ends. 

EDU3> indicates that the Command Window is in normal mode, in MATLAB Student Ver- 
sion. 
^> indicates that the Command Window is in normal mode. 



Command Window 



■+■ n ? x 



<3) New to MATtAB? Watch this Video, see Demos, or read Getting Started . 






This is a Classroom License for instructional use only. 
Research and comir.ercial U3e i3 prohibited. 

MATLAB desktop keyboard shortcuts, such as Ctrl+S, are now customizable . 
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restore previous default settings by following the steps outlined in Hel 
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Figure 1.2: The Command Window. 



1.1.3.2 Command History 

The Command History is a log of the commands we have executed in the command window. 
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Figure 1.3: The Command History. 



1.1.3.3 Workspace 

The workspace consists of a set of variables stored in memory during a MATLAB session. To open 
the Workspace browser, select Desktop > Workspace in the MATLAB desktop, or type 

^> workspace 

at the Command Window prompt. 
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Figure 1.4: Workspace. 



1.1.3.4 Current Folder 

The Current Folder is like the Finder in Mac OS X or Windows Explorer in Windows operating 
systems and allows us to browse through the files and folders. The Current Folder also displays 
details about files in your current directory and within the hierarchy of the folders it contains. 



Current Folder: 



H:\MATLAB 



TEU & 



Figure 1.5: Current Folder. 
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Figure 1.6: Current Folder docked on the desktop. 



1.1.3.5 Start Button 

The MATLAB Start button is located at the lower left corner of the MATLAB desktop and provides 
and easy access to tools, demos, and documentation for the MATLAB installation. 

I 



* Start Ready 



Figure 1.7: Start Button. 
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1.1.3.6 Menu Bar 

The menu bar contains commands for creating, opening, printing, editing, viewing, and manipu- 
lating desktop items. 



I MATLAB 7.11.0 (R2010b) 



File Edit Debug Desktop Window Help 



Figure 1.8: Menu Bar. 



1.1.3.7 Toolbar 

The MATLAB toolbar provides on-screen buttons to access frequently used features such as, copy, 
paste, undo and redo. 
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Figure 1.9: Toolbar. 



1.1.3.8 Keyboard shortcuts 

MATLAB provides keyboard shortcuts for viewing a history of commands and listing contextual 
help. 

1 . The up arrow key, 

2. The tab key, 

3. The semicolon symbol. 



1.1.3.8.1 The Up Arrow Key 

Suppose we want to enter the following equation: 

> y=sin(45) 
But we mistakenly entered 
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> y=sine(45) 

MATLAB returns the following prompt: 

??? Undefined function or method 'sine' for input arguments of type 'double' 

Instead of retyping the equation, press the up arrow key, the mistakenly entered line is displayed. 
Using the left arrow key, move the cursor to the misspelled letter. Make the correction and press 
Return or Enter to execute the command. 

Pressing the up arrow key repeatedly recalls the previously entered commands. Likewise, typ- 
ing the first characters of previously entered line and pressing the up arrow key displays the full 
command line. To execute that line, simply press the Return or Enter key. 

1.1.3.8.2 The Tab Key 

Suppose you forgot how to enter the square root command. Begin typing y=sq in the command 
prompt: 

> y=sq 

Then press the tab key and scroll down to sqrt. Select it and press Return or Enter key. 

> y=sqrt 

1.1.3.8.3 The Semicolon Symbol 

The semicolon symbol at the end of a line suppresses the screen output. This is useful when you 
want to keep your command window clean. 

Type the following entry and press the Return key: 

> y=2+2 

The following output is displayed: 

y = 

4 
Now, press the up arrow key to recall our initial entry 

> y=2+2 

And insert a semicolon as follows: 

> y=2+2; 

No numerical result is displayed however MATLAB stores the value of y in the memory. We can 
recall the value y by simply typing y and pressing Return. 
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1.1.4 MATLAB Help 

MATLAB comes with three forms of online help: help, doc and demos. 

1.1.4.1 Help 

Typing help in the Command Window lists all primary help topics. You can display a topic by 
clicking on the link. 

> help 



muMmtmw 
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User Guides 
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Provides instructions for using help functions, the Help browser, and other resources 

Examples in Documentation 

Lists major examples in the MATLAB documentation 

Programming Tips 

Provides helpful techniques and shortcuts for programming in MATLAB 
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Figure 1.10: Help. 



Or if you know the command or function you need help with, you can type help followed by the 
command or function. For example to learn about clc command, type help clc at the command 
prompt: 



^> help clc 
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Figure 1.11: The output of S> help clc command. 



Also try the following command: ^> help clear 
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MATLAB 7.11.0 (R2D10b) 
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Figure 1.12: The output of S> help clear command. 



To learn about sine function, type help sin at the command prompt: 
S> help sin 

1.1.4.2 Doc 

Obviously, to use help effectively, you need to know what you are looking for. Often times, espe- 
cially when you first start learning an application, it is usually difficult to ask the right questions. 
In the case of MATLAB, doc command is generally better than help. If you type doc in the 
command prompt, MATLAB opens a browser from where you can obtain help easier: 



> doc 
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Figure 1.13: Built-in MATLAB Documentation. 



Like using help sin, try typing doc sin in the command prompt: 
» doc sin 

1.1.4.3 Demos 

You can learn more about MATLAB through demos by typing demo in the command prompt, a list 
of links to demos will open in Help Browser. Demos and online seminars are available at product 
demos and online seminars 9 . 

S> demo 



9 http://www.mathworks.com/products/matlab/demos.html 
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Importing Data from Files [6 min, 36 sec) 



OB Video 



jj 



Working in The Development Environment (4 min, 7 sec) gg Video 



Writing a MATLAB Program (5 min, 45 sec) 



S Video 



zill 



* Start | Ready 



Figure 1.14: Built-in MATLAB Demos. 



1.1.5 Useful Commands and Functions 



For a detailed explanation and examples for each of the following type 'help function' (without 
quotes) at the MATLAB prompt. 
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Command/Function 


Meaning 


clc 


Clear Command Window 


clear 


Remove items from workspace 


who, whos 


List variables in workspace 


workspace 


Display Workspace browser 


cd 


Change working directory 


pwd 


Display current directory 


computer 


Identify information about computer on which MATLAB is running 


ver 


Display version information for Math Works products 


quit 


Terminate MATLAB 


exit 


Terminate MATLAB (same as quit) 



Table 1.1: Useful commands and functions 



1.1.6 Summary of Key Points 

1. MATLAB is a popular technical computing application and Math Works offers a trial version 
of MATLAB on their website, 

2. The MATLAB Desktop consists of Command Window, Command History, Workspace, Cur- 
rent Folder and Start Button, 

3. The up/down arrow keys, the tab key and the semicolon are convenient tools to use the 
Command Window, 

4. MATLAB features an online help, doc and demo, 

5. Various commands and functions make MATLAB experience easier, for example, clc, 
clear and exit. 



1.2 Problem Set 
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Exercise 1.2.1 

Learn about the following terms using help command: 



(Solution on p. 21.) 



1. workspace 

2. plot 

3. clear 



10' 



This content is available online at <http://cnx.org/content/m41463/L2/>. 
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4. format 

5. roots 

Exercise 1.2.2 (Solution on p. 22.) 

List the items found in START button. 

Exercise 1.2.3 (Solution on p. 22.) 

List the items found under DESKTOP menu. 

Exercise 1.2.4 (Solution on p. 23.) 

List the items found under HELP menu. 

Exercise 1.2.5 (Solution on p. 24.) 

Use Function Browser to learn about natural logarithm, (hint: Help Menu > Function 
Browser > Mathematics > Elementary Math > Exponential) 



21 

Solutions to Exercises in Chapter 1 

Solution to Exercise 1.2.1 (p. 19) 

1. 



» help workspace 
WORKSPACE Open Workspace browser to manage workspace 

WORKSPACE Opens the Workspace browser with a view of the variables 
in the current Workspace. Displayed variables may be viewed, 
manipulated, saved, and cleared. 

See also whos, openvar, save. 

Reference page in Help browser 
doc workspace 

> 



» help plot 
PLOT Linear plot. 

PL0T(X,Y) plots vector Y versus vector X. If X or Y is a matrix, 
then the vector is plotted versus the rows or columns of the matrix, 
whichever line up. If X is a scalar and Y is a vector, disconnected 
line objects are created and plotted as discrete points vertically at 
X 



~> help clear 
CLEAR Clear variables and functions from memory. 
CLEAR removes all variables from the workspace . 
CLEAR VARIABLES does the same thing. 
CLEAR GLOBAL removes all global variables. 
CLEAR FUNCTIONS removes all compiled M- and MEX-f unctions . 

CLEAR ALL removes all variables, globals, functions and MEX links. 
CLEAR ALL at the command prompt also removes the Java packages import 
list. 
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~> help format 
FORMAT Set output format. 

FORMAT with no inputs sets the output format to the default appropriate 
for the class of the variable. For float variables, the default is 

FORMAT SHORT. 



» help roots 
ROOTS Find polynomial roots. 

ROOTS (C) computes the roots of the polynomial whose coefficients 
are the elements of the vector C. If C has N+l components, 
the polynomial is C(1)*X~N + ... + C(N)*X + C(N+1). 



Solution to Exercise 1.2.2 (p. 20) 

Following figure illustrates the item found in START button: 



4- 


MATLAB ► 


* 


Toolboxes ► 


* 


Simulink ► 


* 


Blocksets ► 


m 


Shortcuts ► 


& 


Desktop Tools ► 


« 


Web ► 


% 


Get Product Trials 


% 


Check for Updates 


% 


Preferences.,. 


m 


Find Files,., 


% 


Help 


•i 


Demos 


'Start 



Figure 1.15: Start Button 



Solution to Exercise 1.2.3 (p. 20) 

Following figure illustrates the item found under DESKTOP menu: 
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Figure 1.16: Desktop menu items 



Solution to Exercise 1.2.4 (p. 20) 

Following figure illustrates the item found under HELP menu: 
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J MATLAB 7.11.0 (R2010b) 
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Figure 1.17: Help menu items 



Solution to Exercise 1.2.5 (p. 20) 

Following figure shows the solution: 
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I MATLAB 7.: 



File Edit Debug Desktop Window Help 
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B Desktop Tools and Development Environment 
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£3 Rounding and Remainder 

(rS nkrretp Msth 
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Workspace 



<h 



log 

Natural logarithm 



"77 



More Help,. 



The log function operates element-wise on arrays, 
domain includes complex and negative numbers, 
which may lead to unexpected results if used 
unintentionally. 



Y = log(X) returns the natural logarithm of the 
elements of X. For complex or negative Z, where 
2 = x + y l the complex logarithm is returned, 
log ( z) = log (abs ( z) ) + i*atan2 (y, x) 



Figure 1.18: Information about natural logarithm displayed with Search for Functions. 
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Chapter 2 
Getting Started 



2.1 Essentials 




Learning a new skill, especially a computer program in this case, can be overwhelming. However, 
if we build on what we already know, the process can be handled rather effectively. In the preceding 
chapter we learned about MATLAB Graphical User Interface (GUI) and how to get help. Knowing 



This content is available online at <http://cnx.org/content/m41409/Ll/>. 
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the GUI, we will use basic math skills in MATLAB to solve linear equations and find roots of 
polynomials in this chapter. 

2.1.1 Basic Computation 

2.1.1.1 Mathematical Operators 

The evaluation of expressions is accomplished with arithmetic operators as we use them in scien- 
tific calculators. Note the addtional operators shown in the table below: 



Operator 


Name 


Description 


+ 


Plus 


Addition 


- 


Minus 


Subtraction 


* 


Asterisk 


Multiplication 


/ 


Forward Slash 


Division 


\ 


Back Slash 


Left Matrix Division 


A 


Caret 


Power 


* 


Dot Asterisk 


Array multiplication (element-wise) 


./ 


Dot Slash 


Right array divide (element-wise) 


A 


Dot Back Slash 


Left array divide (element-wise) 


A 


Dot Caret 


Array power (element-wise) 



Table 2.1: Operators 



NOTE: The backslash operator is used to solve linear systems of equations, see Sec- 
tion 2.1.5 (Linear Equations). 



important: Matrix is a rectangular array of numbers and formed by rows and columns. 
/ 1 2 3 4 \ 



For example A — 



5 6 7 8 
9 10 11 12 



. In this example A consists of 4 rows and 4 



\ 13 14 15 16 J 

columns and therefore is a 4x4 matrix, (see Wikipedia 2 ). 



http://en.wikipedia.org/wiki/Matrix_%28mathematics%29 
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important: Row vector is a special matrix that contains only one row. In other words, 
a row vector is a lxn matrix where n is the number of elements in the row vector. B 
12 3 4 5 



important: Column vector is also a special matrix. As the term implies, it contains only 
one column. A column vector is an nxl matrix where n is the number of elements in the 
/l\ 



column vector. C = 



2 
3 
4 

\ 5 / 



note: Array operations refer to element-wise calculations on the arrays, for example if x 
is an a by b matrix and y is a c by d matrix then x.*y can be performed only if a=c and 
b=d. Consider the following example, x consists of 2 rows and 3 columns and therefore it 
is a 2x3 matrix. Likewise, y has 2 rows and 3 columns and an array operation is possible. 

1 2 3 \ / 10 20 30 \ / 10 40 90 

a = i I and y — \ ) then x. * y — 



4 5 6 



40 50 60 



160 250 360 



Example 2.1 

The following figure illustrates a typical calculation in the Command Window. 
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MATLAB 7.11.0 (R2D10b) 
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Figure 2.1: Basic arithmetic in the command window. 



2.1.1.2 Operator Precedence 

MATLAB allows us to build mathematical expressions with any combination of arithmetic op- 
erators. The order of operations are set by precedence levels in which MATLAB evaluates an 
expression from left to right. The precedence rules for MATLAB operators are shown in the list 
below from the highest precedence level to the lowest. 

1 . Parentheses () 

2. Power ( A ) 

3. Multiplication (*), right division (/), left division (\) 

4. Addition (+), subtraction (-) 



2.1.2 Mathematical Functions 

MATLAB has all of the usual mathematical functions found on a scientific calculator including 
square root, logarithm, and sine. 
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important: Typing pi returns the number 3.1416. To find the sine of pi, type in sin(pi) 
and press enter. 



important: The arguments in trigonometric functions are in radians. Multiply degrees 
by pi/180 to get radians. For example, to calculate sin(90), type in sin(90*pi/180) . 



warning: In MATLAB log returns the natural logarithm of the value. To find the In of 
10, type in log(10) and press enter, (ans = 2.3026). 



warning: MATLAB accepts loglO for common (base 10) logarithm. To find the log of 
10, type in loglO(lO) and press enter, (ans = 1). 

Practice the following examples to familiarize yourself with the common mathematical functions. 
Be sure to read the relevant help and doc pages for functions that are not self explanatory. 

Example 2.2 

Calculate the following quantities: 

1 



3 2 -r 

2. 5°- 5 -l 

3. frf 2 ford=2 

MATLAB inputs and outputs are as follows: 

1. j|-j- is entered by typing 2~3/ (3~2-l) (ans = 1) 

2. 5 0,5 — 1 is entered by typing sqrt(5)-l (ans = 1.2361) 

3. f d 2 for d=2 is entered by typing pi/4*2~2 (ans = 3.1416) 

Example 2.3 

Calculate the following exponential and logarithmic quantities: 

1. e 2 

2. ln(5 10 ) 

3. loglO 5 

MATLAB inputs and outputs are as follows: 

1. exp(2) (ans = 7.3891) 

2. log((5~10)) (ans= 16.0944) 

3. Iogl0(10~5) (ans = 5) 
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Example 2.4 

Calculate the following trigonometric quantities: 

1. cos(f) 

2. tan (45) 

3. sin (71;) + cos (45) 

MATLAB inputs and outputs are as follows: 

1. cos (pi/6) (ans = 0.8660) 

2. tan(45*pi/180) (ans = 1.0000) 

3. sin(pi)+cos(45*pi/180) (ans = 0.7071) 



2.1.3 The format Function 

The format function is used to control how the numeric values are displayed in the Command 
Window. The short format is set by default and the numerical results are displayed with 4 digits 
after the decimal point (see the examples above). The long format produces 15 digits after the 
decimal point. 

Example 2.5 

Calculate 6 — tan (?) and display results in short and long formats. 

The short format is set by default: 

> theta=(pi/3) 

theta = 

1.0472 

> 

And the long format is turned on by typing format long: 

> theta=(pi/3) 

theta = 

1.0472 

3> format long 
> theta 
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theta = 

1.047197551196598 

2.1.4 Variables 

In MATLAB, a named value is called a variable. MATLAB comes with several predefined vari- 
ables. For example, the name pi refers to the mathematical quantity it, which is approximately pi 

ans = 3.1416 

warning: MATLAB is case-sensitive, which means it distinguishes between upper- and 
lowercase letters (e.g. data, DATA and DaTa are three different variables). Command and 
function names are also case-sensitive. Please note that when you use the command-line 
help, function names are given in upper-case letters (e.g., CLEAR) only to emphasize 
them. Do not use upper-case letters when running functions and commands. 

2.1.4.1 Declaring Variables 

Variables in MATLAB are generally represented as matrix quantities. Scalars and vectors are 
special cases of matrices having size lxl (scalar), lxn (row vector) or nxl (column vector). 

2.1.4.1.1 Declaration of a Scalar 

The term scalar as used in linear algebra refers to a real number. Assignment of scalars in MAT- 
LAB is easy, type in the variable name followed by = symbol and a number: 

Example 2.6 

a = 1 



Command Window -+■ n ? x 






» a=l 



Figure 2.2: Assignment of a scalar quantity. 
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2.1.4.1.2 Declaration of a Row Vector 

Elements of a row vector are separated with blanks or commas. 

Example 2.7 

Let's type the following at the command prompt: 

b = [12 3 4 5] 



Command Window 



» fc=[l 2 3 4 5] 

b = 

12 3 4 5 






Figure 2.3: Assignment of a row vector quantity. 



We can also use the Variable Editor to assign a row vector. In the menu bar, select File 
> New > Variable. This action will create a variable called unnamed which is displayed 
in the workspace. By clicking on the title unnamed, we can rename it to something more 
descriptive. By double-clicking on the variable, we can open the Variable Editor and type 
in the values into spreadsheet looking table. 
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Figure 2.4: Assignment of a row vector by using the Variable Editor. 



2.1.4.1.3 Declaration of a Column Vector 

Elements of a column vector is ended by a semicolon: 

Example 2.8 

c = [1;2;3;4;5;] 
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X 




» c=[l;2;3;4;5; j 


^J 




c = 








1 
2 
3 
4 

5 







Figure 2.5: Assignment of a column vector quantity. 

Or by transposing a row vector with the ' operator: 
c = [12 3 4 5]' 



Command Window 

» c=[l 2 3 4 5] ■ 



-+■ □ ? X 



II 



Figure 2.6: Assignment of a column vector quantity by transposing a row vector with the ' 
operator. 



Or by using the Variable Editor: 
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Figure 2.7: Assignment of a column vector quantity by using the Variable Editor. 



2.1.4.1.4 Declaration of a Matrix 



Matrices are typed in rows first and separated by semicolons to create columns. Consider the 
examples below: 

Example 2.9 

Let us type in a 2x5 matrix: 

d = [2 468 10; 1357 9] 
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Figure 2.8: Assignment of a 2x5 matrix. 
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Figure 2.9: Assignment of a matrix by using the Variable Editor. 



Example 2.10 

This example is a 5x2 matrix: 
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Figure 2.10: Assignment of a 5x2 matrix. 



2.1.5 Linear Equations 

Systems of linear equations are very important in engineering studies. In the course of solving 
a problem, we often reduce the problem to simultaneous equations from which the results are 
obtained. As you learned earlier, MATLAB stands for Matrix Laboratory and has features to 
handle matrices. Using the coefficients of simultaneous linear equations, a matrix can be formed 
to solve a set of simultaneous equations. 

Example 2.11 

Let's solve the following simultaneous equations: 



x + y — 1 



(2.1) 



2x-5y = 9 (2.2) 

First, we will create a matrix for the left-hand side of the equation using the coefficients, 
namely 1 and 1 for the first and 2 and -5 for the second. The matrix looks like this: 



1 1 

2 -5 
The above matrix can be entered in the command window by typing A=[l 1; 2 -5]. 



(2.3) 



Second, we create a column vector to represent the right-hand side of the equation as 
follows: 




(2.4) 
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The above column vector can be entered in the command window by typing B= [1 ; 9] . 

To solve the simultaneous equation, we will use left division operator and issue the fol- 
lowing command: C=A\B. These three steps are illustrated below: 

> A=[l 1; 2 -5] 
A = 

1 1 

2 -5 

> B= [1;9] 

B = 

1 
9 

> C=A\B 
C = 

2 

-1 

> 

The result C indicating 2 and 1 are the values for x and y, respectively. 



2.1.6 Polynomials 

In the preceding section, we briefly learned about how to use MATLAB to solve linear equations. 
Equally important in engineering problem solving is the application of polynomials. Polynomials 
are functions that are built by simply adding together (or subtracting) some power functions, (see 
Wikipedia 3 ). 

ax 2 + bx + c = (2.5) 

f(x) = ax 2 + bx + c (2.6) 



3 http://en. wikipedia.org/wiki/Polynomial 
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The coeffcients of a polynominal are entered as a row vector beginning with the highest power 
and including the ones that are equal to 0. 

Example 2.12 

Create a row vector for the following function: y — 2x 4 + 3x 3 + 5x 2 + x + 10 

Notice that in this example we have 5 terms in the function and therefore the row vector 
will contain 5 elements. p= [2 3 5 1 10] 

Example 2.13 

Create a row vector for the following function: y = 3x 4 + 4x 2 — 5 

In this example, coefficients for the terms involving power of 3 and 1 are 0. The row vector 
still contains 5 elements as in the previous example but this time we will enter two zeros 
for the coefficients with power of 3 and 1: p= [3 4 -5] . 



2.1.6.1 The polyval Function 

We can evaluate a polynomial p for a given value of x using the syntax polyval (p,x) where p 
contains the coefficients of polynomial and x is the given number. 

Example 2.14 

Evaluate f(x) at 5. 

f(x) = 3x 2 + 2x+l (2.7) 

The row vector representing f(x) above is p= [3 2 1] . To evaluate f(x) at 5, we type in: 
polyval (p, 5) . The following shows the Command Window output: 

> p=[3 2 1] 



> polyval (p, 5) 



ans = 

86 
> 
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2.1.6.2 The roots Function 
Consider the following equation: 

ax 2 + bx + c = (2.8) 

Probably you have solved this type of equations numerous times. In MATLAB, we can use the 
roots function to find the roots very easily. 

Example 2.15 

Find the roots for the following: 

0.6x 2 + 0.3* -0.9 = (2.9) 

To find the roots, first we enter the coefficients of polynomial in to a row vector p with 
p= [0 . 6 0.3 -0.9] and issue the r=roots (p) command. The following shows the com- 
mand window output: 

> p=[0.6 0.3 -0.9] 

P = 

0.6000 0.3000 -0.9000 

» r=roots(p) 



r = 



■1.5000 
1.0000 



> 



2.1.7 Splitting a Statement 

You will soon find out that typing long statements in the Command Window or in the the Text 
Editor makes it very hard to read and maintain your code. To split a long statement over multiple 
lines simply enter three periods "..." at the end of the line and carry on with your statement on the 
next line. 

Example 2.16 

The following command window output illustrates the use of three periods: 



ans 
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> sin(pi)+cos(45*pi/180) -sin (pi/2) +cos(45*pi/180)+tan(pi/3) 



2.1463 



> sin(pi)+cos(45*pi/180)-sin(pi/2) 
+cos(45*pi/180)+tan(pi/3) 



ans = 

2.1463 
> 



2.1.8 Comments 

Comments are used to make scripts more "readable". The percent symbol % separates the com- 
ments from the code. Examine the following examples: 

Example 2.17 

The long statements are split to make it easier to read. However, despite the use of de- 
scriptive variable names, it is hard to understand what this script does, see the following 
Command Window output: 

t_water=80; 
t_outside=15; 
inner_dia=0.05; 
thickness=0.006: 
Lambda_steel=48 : 
Alfalnside=2800: 
Alfa0utside=17; 
thickness_insulation=0 . 012 ; 
Lambda_insulation=0 . 03 ; 

r_i=inner_dia/2 

r_o=r_i+thickness 

r_i_insulation=r_o 

r_o_insulation=r_i_insulation+thickness_insulation 

Arealnside=2*pi*r_i 

Area0utside=2*pi*r_o 

Area0utside_insulated=2*pi*r_o_insulation 
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AreaM_pipe=(2*pi*(r_o-r_i))/log(r_o/r_i) 
AreaM_insulation=(2*pi*(r_o_insulation-r_i_insulation)) . . . 

/log(r_o_insulation/r_i_insulation) 
TotalResistance=(l/(AlfaInside*AreaInside))+ . . . 

(thickness/ (Lambda_steel*AreaM_pipe) )+ (1/ (Alf aOutside*AreaOutside) ) 
TotalResistance_insulated=(l/(AlfaInside*AreaInside))+ . . . 

(thickness/ (Lambda_steel*AreaM_pipe))+(thickness_insulation . . . 

/ (Lambda_insulation*AreaM_insulation) ) + ( 1/ (Alf aOutside*AreaOutside_insulated) ) 
Q_dot= (t_water-t_outside) / (TotalResistance*1000) 

Q_dot_insulated=(t_water-t_outside)/(TotalResistance_insulated*1000) 
PercentageReducttion= ( (Q_dot-Q_dot_insulated) /Q_dot) *100 

Example 2.18 

The following is an edited version of the above including numerous comments: 

*/. Problem 16.06 
% Problem Statement 

% Calculate the percentage reduction in heat loss when a layer of hair felt 
% is wrapped around the outside surface (see problem 16.05) 

format short 

% Input Values 

t_water=80; % Water temperature [C] 

t_outside=15; % Atmospheric temperature [C] 

inner_dia=0.05; % Inner diameter [m] 

thickness=0.006; °/„ [m] 

Lambda_steel=48; % Thermal conductivity of steel [W/mK] 

Alfalnside=2800; % Heat transfer coefficient of inside [W/m2K] 

Alfa0utside=17; % Heat transfer coefficient of outside [W/m2K] 

% Neglect radiation 

% Additional layer 

thickness_insulation=0.012; °/ [m] 

Lambda_insulation=0.03; % Thermal conductivity of insulation [W/mK] 



% Output Values 

% Q_dot=(t_water-t_outside)/TotalResistance 

% TotalResistance=(l/(AlfaInside*AreaInside))+(thickness/(Lambda_steel*AreaM))+ . . 

(l/(AlfaOutside*AreaOutside) 

% Calculating the unknown terms 

r_i=inner_dia/2 % Inner radius of pipe [m] 

r_o=r_i+thickness % Outer radius of pipe [m] 
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r_i_insulation=r_o % Inner radius of insulation [m] 

r_o_insulation=r_i_insulation+thickness_insulation °/„ Outer radius of pipe [m] 

Arealnside=2*pi*r_i 

Area0utside=2*pi*r_o 

Area0utside_insulated=2*pi*r_o_insulation 

AreaM_pipe=(2*pi*(r_o-r_i))/log(r_o/r_i) °/ Logarithmic mean area for pipe 

AreaM_insulation=(2*pi*(r_o_insulation-r_i_insulation)) . . . 

/log(r_o_insulation/r_i_insulation) % Logarithmic mean area for insulation 
TotalResistance=(l/(AlfaInside*AreaInside))+ (thickness/ . . . 

(Lambda_steel*AreaM_pipe))+(l/(AlfaOutside*AreaOutside)) 
TotalResistance_insulated=(l/(AlfaInside*AreaInside))+ (thickness/ . . . 

(Lambda_steel*AreaM_pipe))+(thickness_insulation/(Lambda_insulation*AreaM_insulation 

+ (1/ (Alf aOutside*AreaOutside_insulated) ) 
Q_dot=(t_water-t_outside)/(TotalResistance*1000) % converting into kW 

Q_dot_insulated=(t_water-t_outside)/(TotalResistance_insulated*1000) % converting into k 
PercentageReducttion= ( (Q_dot-Q_dot_insulated) /Q_dot) *100 



2.1.9 Basic Operations 



Command 


Meaning 


sum 


Sum of array elements 


prod 


Product of array elements 


sqrt 


Square root 


log 10 


Common logarithm (base 10) 


log 


Natural logarithm 


max 


Maximum elements of array 


min 


Minimum elements of array 


mean 


Average or mean value of arrays 


std 


Standard deviation 



Table 2.2: Basic operations. 
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2.1.10 Special Characters 



Character 


Meaning 


= 


Assignment 


() 


Prioritize operations 


[] 


Construct array 




Specify range of array elements 


? 


Row element separator in an array 


5 


Column element separator in an array 




Continue statement to next line 




Decimal point, or structure field separator 


% 


Insert comment line into code 



Table 2.3: Special Characters 



2.1.11 Summary of Key Points 

1 . M ATLAB has the common functions found on a scientific calculator and can be operated in 
a similar way, 

2. MATLAB can store values in variables. Variables are case sensitive and some variables are 
reserved by MATLAB (e.g. pi stores 3.1416), 

3. Variable Editor can be used to enter or manipulate matrices, 

4. The coefficients of simultaneous linear equations and polynomials are used to form a row 
vector. MATLAB then can be used to solve the equations, 

5. The format function is used to control the number of digits displayed, 

6. Three periods "..." at the end of the line is used to split a long statement over multiple lines, 

7. The percent symbol % separates the comments from the code, anything following % symbol 
is ignored by MATLAB. 



2.2 Problem Set 4 



Determine the value of each of the following. 

Exercise 2.2.1 

6 x 7 + 4 2 - 2 4 



(Solution on p. 49.) 



This content is available online at <http://cnx.org/content/m4 1464/1. 5/>. 
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Exercise 2.2.2 

3 2 +2 3 , 64°' 5 -5 2 



(Solution on p. 49.) 



4 5_ 5 4 -T 4 5 +5 6 +7 8 

Exercise 2.2.3 

^2 i i n5 



sin (270)+ cos (f) 



(Solution on p. 49.) 
Iogl0 z +10 D 

Exercise 2.2.4 
e 2 + 2 3 - In (e 2 ) 

Exercise 2.2.5 
sin(27T)+cos(f) 

Exercise 2.2.6 

tan (§)+ cos (270) 

Exercise 2.2.7 
Solve the following system of equations: 

2x + Ay = 1 
x + 5y = 2 

Exercise 2.2.8 

Evaluate y at 5. 

y — 4x 4 + 3x 2 — x 

Exercise 2.2.9 (Solution on p. 50.) 

Given below is Load-Gage Length data for a type 304 stainless steel that underwent a 
tensile test. Original specimen diameter is 12.7 mm. 5 



(Solution on p. 49.) 



(Solution on p. 49.) 



(Solution on p. 49.) 



(Solution on p. 49.) 



(Solution on p. 49.) 



5 Introduction to Materials Science for Engineers by J. F. Shackelford, Macmillan Publishing Company. ©1985, 
(p.304) 
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Load [kN] 


Gage Length [mm] 


0.000 


50.8000 


4.890 


50.8102 


9.779 


50.8203 


14.670 


50.8305 


19.560 


50.8406 


24.450 


50.8508 


27.620 


50.8610 


29.390 


50.8711 


32.680 


50.9016 


33.950 


50.9270 


34.580 


50.9524 


35.220 


50.9778 


35.720 


51.0032 


40.540 


51.816 


48.390 


53.340 


59.030 


55.880 


65.870 


58.420 


69.420 


60.960 


69.670 (maximum) 


61.468 


68.150 


63.500 


60.810 (fracture) 


66.040 (after fracture) 



Table 2.4 



The engineering stress is defined as a — j, where P is the load [N] on the sample with an 
original cross-sectional area A [m 2 ] and the engineering strain is defined as e — y, where 
Al is the change in length and / is the initial length. 



Compute the stress and strain values for each of the measurements obtained in the 
tensile test. 
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Solutions to Exercises in Chapter 2 

Solution to Exercise 2.2.1 (p. 46) 

> (6*7) +4~2-2~4 (ans = 42) 
Solution to Exercise 2.2.2 (p. 46) 

> ( (3~2+2~3) / (4~5-5~4) ) + ( (sqrt (64) -5~2) / (4~5+5~6+7~8) ) (ans = 0.0426) 
Solution to Exercise 2.2.3 (p. 47) 

> logl0(10~2)+10~5 (ans = 100002) 
Solution to Exercise 2.2.4 (p. 47) 

> exp(2)+2~3-log(exp(2)) (ans = 13.3891) 
Solution to Exercise 2.2.5 (p. 47) 

> sin(2*pi)+cos(pi/4) (ans = 0.7071) 
Solution to Exercise 2.2.6 (p. 47) 

> tan(pi/3)+cos(270*pi/180)+sin(270*pi/180)+cos(pi/3) (ans = 1.2321) 
Solution to Exercise 2.2.7 (p. 47) 



> A=[2 4; 1 5] 

A = 

2 4 

1 5 

> B=[l; 2] 

B = 

1 
2 

> Solution=A\B 

Solution = 

-0.5000 
0.5000 



Solution to Exercise 2.2.8 (p. 47) 



> p=[4 3 -1 0] 
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> polyval(p,5) 



ans = 



2570 



> 



Solution to Exercise 2.2.9 (p. 47) 

First, we need to enter the data sets. Because it is rather a large table, using Variable Editor is 
more convenient. See the figures below: 



|TFf[ Variable Editor - Load_M ■+■ D ? X, 


Workspace 


■+■ n ? 


X 


3g 


a 


si 


f. • 






E^? No valid pi... •* 


D ^J r X 


J ^j "jd 9 1 1 E^J> Select data to .. . » 




ffl 


Load_M <21>1 doubled 


Name L. 


Value 


Hi 




1 


2 


3 


4 


5 


6 




3]l.oad_N <Zlxl doubled 





i 















* 


2 


4890 














3 


9779 














4 


14670 














5 


19560 














6 


24450 














7 


27620 












— ' 


S 


29390 














9 


32680 














10 


33950 














11 


34580 














12 


35220 














<l 1 


±J 


13 


35720 


















14 


40540 














Command History "+ 1 D ? 


X 


15 


48390 














* — 30/10/2011 12:03 — * 




16 


590.30 














17 


65870 












'1 


<i i >r 



Figure 2.11: Load in Newtons 
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r^l Variable Editor - Length_mm ■+■ □ ? 


X 


Workspace ■+■ □ * X 


*k 


i 


& 


d ' 


% 


J ^ Select d... » 


□ zl * 


X 


a] |-6j Qb] ^ Bjjj II [~/TJ p|ot(Lengtri_mm) t 


a 


Lengthjnm <21xl double> 


Name £. 


Value Mi 




1 


2 


3 


4 


5 


6 


±1 


J] Lengthjnm <21xl double> 5C 
jLoadJI <21xldouble> 


6 


50.8508 












7 


50.8610 












3 


50.8711 












9 


50.9016 












10 


50.9270 












11 


50.9524 












12 


50.9778 












13 


51.0032 












H 


51.8160 












15 


53.3400 












16 


55.8800 












17 


58.4200 












<l 1 >l 


18 


60.9600 














19 


61.4680 












Command History "+ 1 D ^ X 


20 


63.5000 












% — 30/10/2011 12:03 — % 


21 


66.0400 




























<l 1 







Figure 2.12: Extension length in mm. 



Next, we will calculate the cross-sectional area. 

Area=pi/4* (0.0127~2) 

Area = 

1.2668e-004 
Now, we can find the Stress values with the following, note that we are obtaining results in MPa: 
Sigma= (Load.N . /Area) *1(T (-6) 

Sigma = 




38.6022 
77.1964 
115.8065 
154.4086 
193.0108 
218.0351 
232.0076 
257.9792 
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268.0047 
272.9780 
278.0302 
281.9773 
320.0269 
381.9955 
465.9888 
519.9844 
548.0085 
549.9820 
537.9830 
480.0403 

For strain calculation, we will first find the change in length: 

Delta_L=Length_mm-50 . 800 
Delta.L = 



0.0102 

0.0203 

0.0305 

0.0406 

0.0508 

0.0610 

0.0711 

0.1016 

0.1270 

0.1524 

0.1778 

0.2032 

1.0160 

2.5400 

5.0800 

7.6200 
10.1600 
10.6680 
12.7000 
15.2400 

Now we can determine Strain with the following: 
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Epsilon=Delta_L . /50 . 800 

Epsilon = 


0.0002 
0.0004 
0.0006 
0.0008 
0.0010 
0.0012 
0.0014 
0.0020 
0.0025 
0.0030 
0.0035 
0.0040 
0.0200 
0.0500 
0.1000 
0.1500 
0.2000 
0.2100 
0.2500 
0.3000 

The final results can be tabulated as foolows: 

[Sigma Epsilon] 
ans = 













38 


,6022 


0. 


,0002 


77 


,1964 


0. 


.0004 


115. 


,8065 


0. 


,0006 


154. 


,4086 


0. 


,0008 


193. 


,0108 


0. 


,0010 


218. 


,0351 


0. 


,0012 


232 


,0076 





,0014 


257. 


,9792 


0. 


,0020 
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268.0047 0.0025 

272.9780 0.0030 

278.0302 0.0035 

281.9773 0.0040 

320.0269 0.0200 

381.9955 0.0500 

465.9888 0.1000 

519.9844 0.1500 

548.0085 0.2000 

549.9820 0.2100 

537.9830 0.2500 

480.0403 0.3000 



Chapter 3 
Graphics 



3.1 Plotting in MATLAB 




A picture is worth a thousand words, particularly visual representation of data in engineering is 
very useful. MATLAB has powerful graphics tools and there is a very helpful section devoted to 
graphics in MATLAB Help: Graphics. Students are encouraged to study that section; what follows 



This content is available online at <http://cnx.org/content/m4 1442/1. 2/>. 
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is a brief summary of the main plotting features. 

3.1.1 Two-Dimensional Plots 

3.1.1.1 The plot Statement 

Probably the most common method for creating a plot is by issuing plot (x , y) statement where 
function y is plotted against x. 

Example 3.1 

Type in the following statement at the MATLAB prompt: 

x=[-pi: . l:pi] ; y=sin(x); plot(x,y); 
After we executed the statement above, a plot named Figure 1 is generated: 



0.8 


I I I I I --- V I I 


0.6 


/ \ 


0.4 


/ \ 


0.2 


/ \ 





- X / 


0.2 


\ / 


0.4 


\ / 


0.6 


\ / 


0.8 
.1 


i i X__-^ i i i i i 



-4 



-2 



Figure 3.1: Graph of sin(x) 
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Having variables assigned in the Workspace, x and y=sin(x) in our case, we can also select x and 
y, and right click on the selected variables. This opens a menu from which we choose plot(x,y). 
See the figure below. 



1 Workspace -" □ ? x 


i§ E3 *| « Bi m plotfey) T 


Name s - 


Value 1 


EEH < 1x63 doubles 




Open Selection Ctrl+D 




Save As.,. 




Copy Ctrl+C 


Duplicate 




Delete Delete 


1\ 


Rename 


= Edit Value 




IB plotfay) 


"J 


Plot as multiple series 


'r 


barfoy) 




area(x,y) 




r= piefay) 




More Plots.,. 




Plot Catalog... 


hplot (x, y, ■ DisplayNaite ■ , ■ y 




L clc T | 


^i i >r 



Figure 3.2: Creating a plot from Workspace. 



3.1.1.2 Annotating Plots 

Graphs without labels are incomplete and labeling elements such as plot title, labels for x and 
y axes, and legend should be included. Using up arrow, recall the statement above and add the 
annotation commands as shown below. 

x=[-pi: . l:pi] ;y=sin(x) ;plot(x,y) ; title ('Graph of y=sin(x) ') ;xlabel('x') ;ylabel( 'sin( 

Run the file and compare your result with the first one. 
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1 

0.8 
0.6 
0.4 
0.2 
f 

en 

-0.2 
-0.4 
-0.6 
-0.8 








Graph of 


y=sin(x) 












1 s* "-v 
















































































































X_-/ i 


1 







Figure 3.3: Graph of sin(x) with Labels. 



aside: Type in the following at the MATLAB prompt and learn additional commands to 
annotate plots: 



help gtext 
help legend 
help zlabel 



3.1.1.3 Superimposed Plots 

If you want to merge data from two graphs, rather than create a new graph from scratch, you can 
superimpose the two using a simple trick: 
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'/. This script generates sin(x) and cos(x) plot on the same graph 
"/. initialize variables 

x=[-pi: . l:pi] ; "/.create a row vector from -pi to +pi with .1 increments 
yO=sin(x); "/.calculate sine value for each x 
yl=cos(x); "/calculate cosine value for each x 
% Plot sin(x) and cos(x) on the same graph 
plot(x,yO,x,yl) ; 

title ('Graph of sin(x) and cos(x)'); "/Title of graph 
xlabel('x'); "/Label of x axis 

ylabeK'sin(x) , cos(x)'); "/Label of y axis 

legendCsin(x) ' , 'cos(x) ') ; "/Insert legend in the same order as yO and yl calculated 
grid on "/Graph grid is turned 



Graph ofsin(x) and cos(x) 




Figure 3.4: Graph of sin(x) and cos(x) in the same plot with labels and legend. 
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3.1.1.4 Multiple Plots in a Figure 

Multiple plots in a single figure can be generated with subplot in the Command Window. How- 
ever, this time we will use the built-in Plot Tools. Before we initialize that tool set, let us create the 
necessary variables using the following script: 

"/ This script generates sin(x) and cos(x) variables 

clc °/ Clears command window 

clear all "/.Clears the variable space 

close all °/ Closes all figures 

Xl=[-2*pi: .l:2*pi] ; "/.Creates a row vector from -2*pi to 2*pi with .1 increments 

Yl=sin(Xl) ; "/Calculates sine value for each x 

Y2=cos(Xl); "/Calculates cosine value for each x 

Y3=Y1+Y2; "/Calculates sin(x)+cos(x) 

Y4=Y1-Y2; "/Calculates sin(x)-cos(x) 

Note that the above script clears the command window and variable workspace. It also closes 
any open Figures. After running the script, we will have XI, Yl, Y2, Y3 and Y4 loaded in the 
workspace. Next, select File > New > Figure, a new Figure window will open. Click "Show Plot 
Tools and Dock Figure" on the tool bar. 
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1*/ Figure 1 


- In 


X 


File Edit View Insert Tools Desktop Window Help 


Tl 


Q S a ^ k \ s. O ® * A T S DB no 





X" 



Workspace 



=| E' *l H f 
I.. 



Show Plot Tools and Dock Figure — 
rrrai 



3Y1 
3Y2 
3Y3 
^4 



Command Histi 



tan [pi/3 

clc 

clear 

- i — 04/10/; 

clear 
clc 

X 

clc 



Figure 3.5: Plot Tools 



Under New Subplots > 2D Axes, select four vertical boxes that will create four subplots in one 
figure. Also notice, the five variables we created earlier are listed under Variables. 
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File Edit View Insert Tools Debug Desktop Window Help * | ? x 


D a A ^ \[k\ \ %■ r> $ v -G ^ t | ^ | □ El | n d ffl m b & \n 


Figure Palette ■■■■■■ '*- H ? X 


■ ■ ■ 


.-.■■■ -+' n ? x 


T New Subplots 






L 


G 2D Axes |fflj 
^ 3D Axes EB ► 


■nnn 
■nnn 
■nnn 
■nnn 
nnnn 




▼Variables 


EBXl 1x126 
EBYl 1x126 
EBY2 1x126 
fflY3 1x125 
EBY4 1x126 


Cancel 






T Annotations 


■ 


\ Line 


\ Arrow 






\ Double Arrow 






^ Text Arrow 






T Text Box 






□ Rectangle 






O Ellipse 








Add Data... | 


Property Editor | 
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Figure 3.6: Creating four sub plots. 



After the subplots have been created, select the first supblot and click on "Add Data". In the 
dialog box, set X Data Source to XI and Y Data Source to Yl. Repeat this step for the remaining 
subplots paying attention to Y Data Source (Y2, Y3 and Y4 need to be selected in the subsequent 
steps while XI is always the X Data Source). 



63 



File Edit View Insert Tools Debug Desktop Window Help 
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Figure 3.7: Adding data to axes. 



Next, select the first item in "Plot Browser" and activate the "Property Editor". Fill out the fields 
as shown in the figure below. Repeat this step for all subplots. 
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Figure 3.8: Using "Property Editor". 



Save the figure as sinxcosx . fig in the current directory. 
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Figure 3.9: The four subplots generated with "Plot Tools". 
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) Figures- Figure 1 


_ 


n 


X 


File Edit View Insert Tools Debug Desktop Window Help 


* | r x 


Q £3 J £ fe \. , n ® * ^ - S DS u o 


ffl m b s^~ 



Graph of sin(x) 




-2 2 

x 

Graph of sin(x)+cos(x) 




-2 2 

x 

Graph of sin(x)-cos[x) 




Figure 3.10: The four subplots in a single figure. 



3.1.2 Three-Dimensional Plots 

3D plots can be generated from the Command Window as well as by GUI alternatives. This time, 
we will go back to the Command Window. 
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3.1.2.1 The plot3 Statement 

With the X1,Y1,Y2 and Y2 variables still in the workspace, type in plot3(Xl ,Y1 ,Y2) at the 
MATLAB prompt. A figure will be generated, click "Show Plot Tools and Dock Figure". 
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Figure 3.11: A raw 3D figure is generated with plot3. 



Use the property editor to make the following changes. 
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Figure 3.12: 3D Property Editor. 



The final result should look like this: 
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*/ Figures - Figure 1 
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q a a ^ k % i © ® « & * & □ b 'Q 


ffl m b s^~ 



Graph of x, sin(x), cos(x) 




Figure 3.13: 3D graph of x, sin(x), cos(x) 



Use help or doc commands to learn more about 3D plots, for example, image (x) , surf (x) and 
mesh(x). 

3.1.3 Generate Code 



A code can be generated to reproduce the plots. To initialize this process, recall sinxcosx.fig 
and select File > Generate Code. 
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Figure 3.14: Generating code to reproduce a plot. 
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Figure 3.15: M-Code generation in progress. 



function createf igure2(Xl , Yl, Y2, Y3, Y4) 
•/.CREATEFIGURE2(X1,Y1,Y2,Y3,Y4) 



*/. XI 

I Yl 

*/. Y2 

*/. Y3 

'/. Y4 



vector of x data 

vector of y data 

vector of y data 

vector of y data 

vector of y data 



% Auto-generated by MATLAB on 05-0ct-2011 12:43:49 

% Create figure 
figurel = figure; 

% Create axes 

axesl = axesCParent' ,f igurel , 'YGrid' , 'on' , 'XGrid' , 'on' , . . . 

'Position' ,[0.13 0.791155913978495 0.775 0.11741935483871]); 
box(axesl, 'on') ; 
hold(axesl,'all'); 

% Create title 

title ('Graph of sin(x)'); 

7. Create xlabel 
xlabel('x'); 
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7, Create ylabel 
ylabeK'Sin(x)'); 

% Create plot 

plot (XI ,Y1, 'Parent ', axesl, 'DisplayName' , 'Yl vs XI'); 

% Create axes 

axes2 = axes( 'Parent ' ,figurel, 'YGrid' , 'on' , 'XGrid' , 'on' ,.. . 

'Position' ,[0.13 0.572069892473118 0.775 0.11741935483871]); 
box(axes2, 'on') ; 
hold (axes2, 'all'); 

% Create title 

title ('Graph of cos(x)'); 

'/, Create xlabel 
xlabel('x'); 

'/, Create ylabel 
ylabeK'Cos(x)'); 

% Create plot 

plot (XI, Y2, 'Parent', axes2, 'DisplayName' , 'Y2 vs XI'); 

% Create axes 

axes3 = axesCParent' ,figurel, 'YGrid' , 'on' , 'XGrid' , 'on' ,.. . 

'Position', [0.13 0.352983870967742 0.775 0.11741935483871]); 
box(axes3, 'on') ; 
hold(axes3,'all'); 

% Create title 

titleCGraph of sin(x)+cos(x) ') ; 

'/, Create xlabel 
xlabel('x'); 

'/, Create ylabel 
ylabel('Cos(x)+Sin(x)'); 

% Create plot 

plot (XI, Y3, 'Parent', axes3, 'DisplayName' , 'Y3 vs XI'); 
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% Create axes 

axes4 = axes( 'Parent ' ,figurel, 'YGrid' , 'on' , 'XGrid' , 'on' ,.. . 

'Position' ,[0.13 0.133897849462366 0.775 0.11741935483871]); 
box(axes4, 'on') ; 
hold(axes4, 'all') ; 

% Create title 

titleCGraph of sin(x)-cos(x) ') ; 

% Create xlabel 
xlabel('x'); 

'/, Create ylabel 
ylabel('Sin(x)-Cos(x)'); 

% Create plot 

plot(Xl,Y4, 'Parent' ,axes4, 'DisplayName' , 'Y4 vs XI'); 

As you can see, the file assumes you are using the same variables originally used to create the 
graph, therefore the variables need to be passed as arguments in the future executions of the gen- 
erated code. 

3.1.4 Summary of Key Points 

1. plot (x, y) and plot3 (XI ,Y1 , Y2) statements create 2- and 3-D graphs respectively, 

2. Plots at minimum should contain the following elements: title, xlabel, ylabel and 
legend, 

3. Annotated plots can be easily generated with GUI Plot Tools, 

4. MATLAB can generate code to reproduce plots. 



3.2 Problem Set 2 

Exercise 3.2.1 (Solution on p. 77.) 

Plot y — a + bx, using the specified coefficients and ranges (use increments of 0.1): 

a. a = 2, b = 0.3 for < x < 5 

b. a = 3, b = for < x < 10 

c. a = 4, b = -0.3 for < x < 15 

Exercise 3.2.2 (Solution on p. 77.) 

Plot the following functions, using increments of 0.01 and a — 6, b — 0.8, < x < 5: 



This content is available online at <http://cnx.org/content/m4 1466/1. 6/>. 
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a. y — a + x 

b. y — ax 

c. y — asm(x) 

Exercise 3.2.3 

Plot function y — — 

Exercise 3.2.4 

Data collected from Boyle's Law experiment are as follows: 



for y^K < x < IOtt using increments of 



100 



(Solution on p. 79.) 
(Solution on p. 80.) 



Volume [cm A 3] 


Pressure [Pa] 


7.34 


100330 


7.24 


102200 


7.14 


103930 


7.04 


105270 


6.89 


107400 


6.84 


108470 


6.79 


109400 


6.69 


111140 


6.64 


112200 



Table 3.1 

Plot a graph of Pressure vs Volume, annotate your graph. 

Exercise 3.2.5 

The original data collected from Boyle's 3 experiment are as follows: 



(Solution on p. 81.) 



Volume [tube-inches] 


Pressure [inches-Hg] 


12 


29.125 


10 


35.000 


8 


43.688 


6 


58.250 


5 


70.000 


4 


87.375 


3 


116.500 



3 Introduction to Engineering: Modeling and Problem Solving by J. B. Brockman, John Wiley and Sons, Inc. 
©2009, (p.246) 
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Table 3.2 



(Solution on p. 82.) 



Plot a graph of Pressure vs Volume, annotate your graph. 

Exercise 3.2.6 

Display the two plots created earlier in one plot. 

Exercise 3.2.7 (Solution on p. 83.) 

A tensile test of SAE 1020 steel produced the data below 4 experiment are as fol- 
lows: 



Extension [mm] 


Load [kN] 


0.00 


0.0 


0.09 


1.9 


0.31 


6.1 


0.47 


9.4 


2.13 


11.0 


5.05 


11.7 


10.50 


12.0 


16.50 


11.9 


23.70 


10.7 


27.70 


9.3 


34.50 


8.1 



Table 3.3 



Plot a graph of Load vs Extension, annotate your graph. 

Exercise 3.2.8 (Solution on p. 84.) 

Given below is Stress-Strain data for a type 304 stainless steel. 5 
follows: 



experiment are as 



4 Introduction to Materials Science for Engineers I Instructor's Manual by J. F. Shackelford, Macmillan Publishing 
Company. ©1992, (p.440) 

5 Introduction to Materials Science for Engineers by J. F. Shackelford, Macmillan Publishing Company. ©1985, 
(p.304) 
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Stress [MPa] 


Strain [mm/mm] 


0.0 


0.0000 


38.6 


0.0002 


77.2 


0.0004 


115.8 


0.0006 


154.4 


0.0008 


193.0 


0.0010 


218.0 


0.0012 


232.0 


0.0014 


258.0 


0.0020 


268.0 


0.0025 


273.0 


0.0030 


278.0 


0.0035 


282.0 


0.0040 


320.0 


0.0200 


382.0 


0.0500 


466.0 


0.1000 


520.0 


0.1500 


548.0 


0.2000 


550.0 


0.2100 


538.0 


0.2500 


480.0 


0.3000 



Table 3.4 



Plot a graph of Stress vs Strain, annotate your graph. 
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Solutions to Exercises in Chapter 3 

Solution to Exercise 3.2.1 (p. 73) 

a. 

a=2; b=.3; x=[0:.l:5]; y=a+b*x; 
plot (x , y) , title ('Graph of y=a+bx') ,xlabel( 'x') ,ylabel('y') ,grid 



a=3; b=.0; x=[0:.l:10]; y=a+b*x; 
plot (x , y) , title ('Graph of y=a+bx') ,xlabel( 'x') ,ylabel('y') ,grid 



a=2; b=.3; x=[0:.l:5]; y=a+b*x; 
plot (x , y) , title ('Graph of y=a+bx') ,xlabel( 'x') ,ylabel('y') ,grid 

Solution to Exercise 3.2.2 (p. 73) 

a. 

a=6; b=.8; x=[0:.01:5]; y=a+x."b; 
plot (x , y) , title ( ' Graph of y=a+x~b ' ) , xlabel ( ' x ' ) , ylabel ( ' y ' ) , grid 



10 

9.5 

9 

8.5 

>-. 8 

7.5 

7 

6.5 
c 






Graph ol 


y=a+x b 
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a=6; b=.8; x=[0:.01:5]; y=a*x.~b; 
plot (x , y) , title ('Graph of y=ax~b') ,xlabel( 'x') ,ylabel('y') ,grid 



25 



20 



15 - 



Graph of y=ax b 



10 -■ 






0.5 1 1.5 



2.5 



3.5 



4.5 



a=6; x=[0:.01:5]; y=a*sin(x); 
plot(x,y) , title ('Graph of y=a*sin(x) ') ,xlabel('x') ,ylabel('y') ,grid 
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Graph of j 


f=a*sin(x) 
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Solution to Exercise 3.2.3 (p. 74) 



x = pi/100:pi/100:10*pi; 
y = sin(x) ./x; 
plot(x,y) , title ( 'Graph of y=sin(x)/x') ,xlabel('x') ,ylabel('y') ,grid 
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*/ Figure 1 



File Edit View Insert Tools Desktop Window Help 



n x 



Q a a ^ k ^^^®i^- S DS QO 



□ .8 



□ .6 



□ .4 



□ .2 



-□.2 



-□.4 



Graph of y=sin[x)/x 






1D 



15 



20 
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Figure 3.16: Graph of y 



sin(x) 



Solution to Exercise 3.2.4 (p. 74) 



Pressure= [100330 , 102200 , 103930 , 105270 , 107400 , 108470 , 109400 ,111140 ,112200] ; 
Volume= [7. 34 ,7. 24, 7. 14, 7. 04, 6. 89 ,6. 84 ,6. 79 ,6. 69 ,6. 64] ; 
plot(Volume, Pressure) , title ( 'Pressure Volume Graph') ,xlabel( 'Volume') ,ylabel(' Pressure' 
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1.14 



1.12 



1.1 



1.08 



x 10 



Pressure Volume Graph 



CO 

•Xi 
(D 



°- 1.06 - 



1.04 



1.02 






1 

6.6 6.7 6.8 6.9 7 7.1 7.2 7.3 7.4 

Volume 
Solution to Exercise 3.2.5 (p. 74) 



> P=[29. 125,35,43. 688,58. 25,70,87. 375,116. 5] ; 
> V=[12,10,8,6,5,4,3]; 
;» plot(V,P) ,title( 'Pressure Volume Graph') ,xlabel( 'Volume') ,ylabel( 'Pressure') , grid 



82 



CHAPTER 3. GRAPHICS 



120 
110 
100 
90 
80 
70 
60 
50 
40 
30 
20 



Pressure Volume Graph 






7 8 

Volume 



10 



11 



12 



Solution to Exercise 3.2.6 (p. 75) 
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120 
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Volume 
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Solution to Exercise 3.2.7 (p. 75) 



Extension [0.00 ,0.09 ,0.31 ,0.47 ,2. 13 ,5. 05 ,10. 50 ,16. 50 ,23. 70 ,27. 70 ,34. 50] ; 
Load= [0.0, 1.9, 6. 1,9. 4, 11. 0,11. 7, 12. 0,11. 9, 10. 7, 9. 3, 8.1] ; 



plot (Extension, Load) , title ( 'Load versus Extension Curve') ,xlabel( 'Extension') ,ylabel('L 
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Solution to Exercise 3.2.8 (p. 75) 

The data can be entered using Variable Editor: 



H Variable Editor - Stress ■+■ □ ? x 


ffl-a 


A 


* 
& 


A - 


% 


Base 


E^P No valid pi,. , » 


n z. 


? X 


2 Stress <21x 1 double> 




1 


2 


3 


4 


5 


S 




1 















* 


2 


38.6022 














3 


77. 196-1 














4 


115.8065 














5 


154.4086 














6 


193.0108 














7 


218.0351 














3 


232.0076 












■^^ 


9 


257.9792 





























1 Workspace 



2 Strain 
distress 



Hj [T*j Qb] ajfl B ^3 Select data to 



Name L. 



Value 



<21xldouble> 
<21xl double > 



Then execute the following: 



plot (Strain, Stress) , title ('Stress versus Strain Curve') ,xlabel( 'Strain [mm/mm] ') ,yla 
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Stress versus Strain Curve 
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Chapter 4 

Introductory Programming 



4.1 Writing Scripts to Solve Problems 1 




MATLAB provides scripting and automation tools that can simplify repetitive computational tasks. 
For example, a series of commands executed in a MATLAB session to solve a problem can be saved 
in a script file called an m-file. An m-file can be executed from the command line by typing the 



This content is available online at <http://cnx.org/content/m4 1440/1. 2/>. 
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name of the file or by pressing the run button in the built-in text editor tool bar. 

4.1.1 Script Files 

A script is a file containing a sequence of MATLAB statements. Script files have a filename 
extension of .m. By typing the filename at the command prompt, we can run the script and obtain 
results in the command window. 



I MATLAB 7.11.0 (R2D10b) 



File Edit View Debug Desktop Window Help 

D £3 | H ^ u |J*rf^|t»,%| Current Folder:| H:\MATLAB\M-Files_HeatTransfer 

Shortcuts J How to Add What's New 



D I X 



lDJ a 



Current Folder 


H- □ ? X 


'... « M-File... ► 


* p\ © *- 


JName t- 



__ : Archive 

..; html 

53 Probleml6_01m 

Q Probleml6jH.m 

DProbleml6_Q3,m 

D Probleml6_04.m 

D Probleml6_05,asv 

Q ProblemlSiK.m 

J Probleml6_Q6,asv 

Q ProblemlSiK.m 

D ThermalConductivity.asv 

B ThermalConductivity.m 



alCcnductivrty.rn i;MATLAB S 



Stainless Steel 



Command Window 



isirfcda ?.l'jrc,iri'jnL = 



199.7688 






Workspace 



■i cj *l a ■ 



Name ^~ 



Value 



3 deltaT_12 
3 deltaT_23 
3 deltaT_45 
3 deltaT_56 
3 deltaT_AI_aver... 
3 deltaT_SS_ave... 
3 deltaZJ 
3 deltaZ_2 
3 deltaZ_3 



40 

40 

8.6000 

8.7000 

8.6500 

40 

1 

1 

3 



J 



Command History 



27/09/2011 11:39 



# Start 



Figure 4.1: Number of m-files are displayed in the Current Folder sub-window. 



A sample m-file named ThermalConductivity.m is displayed in Text Editor below. Note the 
triangle (in green) run button in the tool bar, pressing this button executes the script in the command 
window. 
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- h:\V1ATLAB\V-=i es_HeatTransfer\Thf 



Edit Text Go Cell Tools Debug Desktop Window Help 



JLl 



' | ? X 



-■ 



£_ ' | & ^| « * + ^ | S -i £) <B *B HP ii *■ | Steclef^ []| h 



m b en 



B Open file (Ctrl +0) I 
i r' u i 



1.1 



■;■-. 



,% 0» 



% 5tainless Steel 

lambda 5tairilessSteel=14 . 4; 

delfcaZ_l=i; 

deltaZ_2=l; 

deltaZ_SS_average=[deltaZ_l+deltaZ_2>/2; 

t 1=4 40; 

t2=40Q; 

t;3=3€Q; 

deltaT_12=tl-t2; 

deltaT_23=t2-t3; 

deltaT_S5_average=[deltaT_12+deltaT_23J/2; 



% Al'jir.iri'jir. 

deltaZ 3=3; 

deltaZ_4=3; 

deltaZ_Al_average= 

t4=270; 

t5=261.4; 

t6=252.7; 

deltaT_45=t4-t5; 

deltal_56=t5-t6; 

deltaT_Al_average= 



Thermal Conductivity of Staninleaa Steel [W/mK] 

Distance between the first set of thermocouples [cm] 

Distance between the second set of thermocouples [cm] 

Average Distance [era] 

Temperature 1 [C] 

Temperature 2 [C] 

Temperature 3 [C] 

Temperature difference between trjerm.ccouples 

Temperature difference betwe 

Average temperature differ 



and 2 [CI 

r.r.ezrmccouples 2 and 3 [C] 
for Stainless Steel [C] 



ideltaZ 3+deltaZ 4>/2; 



ideltaT 45+deltaT 56) /2; 



% Distance between the third set of thermocouples [cm] 

\ Distance between the fourth set of thermocouples [cm] 

\ Average Distance [era] 

% Temperature 4 [C] 

\ Temperature 5 [C] 

% Temperature 6 [C] 

% Temperature difference between thermocouples 4 and 5 

% Temperature difference between thermocouples 5 and 6 

h Average temperature difference for aluminum [C] 



% Calculation for Thermal Conductivity of Aluminum [W/raK] 
lamfcda_ATjmin-am=lambda_5tainlesaSt:eel* [deltaT_33_ave~ace/deltaT_Al_ave~ace) ' ( de 1 t aZ_&l_ 



age/del taZ_SS_ 



| script 



Ln 1 Col 1 



Figure 4.2: The content of ThermalConductivity .m file is displayed in Text Editor. 



Now let us see how an m-file is created and executed. 

Example 4.1 

A cylindrical acetylene bottle with a radius r=0.3 m has a hemispherical top. The height 
of the cylindrical part is h=1.5 m. Write a simple script to calculate the volume of the 
acetylene bottle. 

To solve this problem, we will first apply the volume of cylinder equation (4.1). Using the 
volume of sphere equation (4.2), we will calculate the volume of hemisphere (4.3). The 
total volume of the acetylene bottle is found with the sum of volumes equation (4.4). 



'cylinder 



7lr 2 h 



(4.1) 



•'sphere — ~ ^ r 



^top = ^ Kr 



(4.2) 
(4.3) 
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''acetylene bottle — ''cylinder T nop (4.4) 

To write the script, we will use the built-in text editor. From the menu bar select File > 
New > Script. The text editor window will open in a separate window. First save this 
file as AcetyleneBottle .m. In that window type the following code paying attention to 
the use of percentage and semicolon symbols to comment out the lines and suppress the 
output, respectively. 

% This script computes the volume of an acetylene bottle with a radius r=0.3 m, 
% a hemispherical top and a height of cylindrical part h=1.5 m. 
r=0.3; 
h=1.5; 

Vol_top= (2*pi*r~3) /3 ; 
Vol_cyl=pi*r~2*h; 
Vol_total=Vol_top+Vol_cyl 



°/ Radius [m] 

*/. Height [m] 

% Calculating the volume of hemispherical top [m3] 

% Calculating the volume of cylindrical bottom [m3] 

*/, Calculating the total volume of acetylene bottle [m3] 



I l«!U..K.FM.!^!ffl«fl!JBIMJIi:mi 



n x 



File Edit Text Go Cell Tools Debug Desktop Window Help 



<|* x 



3 A 



n a 



* 



*J C"|^^^|»4i*^|Sfej'^l^lCl*DflDJD 



Stack: Base 



in b &\n 



- ID 



1.1 



% 



% This script, compute 

% and a height of cyl 

r=G.3; 

h=l . 5„- 

Vol_cop= ( 2 *pi *r"3 \ /3 ; 

Vol_cyl=pi *r*2 *"h; 

Val_tctal=Vol_tap+Val_cyl 



etyle 



r.err.^spr: 






drical part h=l . 5 : 



% Radius [m] 

% rieight [ml 

% Calculating the volume of hemispherical top [m31 

\ Calculating the volume of cylinderical bottom [m3] 

% Calculating the total volume of acetylene battle [m3] 



| script 



Ln 7 Col 1DQ OVR 



Figure 4.3: Script created with the built-in text editor. 
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After running the script by pressing the green button in the Text Editor tool bar, the output 
is displayed in the command window as shown below. 



File Edit View Debug Desktop Window | Help 



D I X 



Q6 «ll^ kE^l « 



Current Folder: H:\MATLAB\M-Files_HeatTransfer 



11 A a 



Shortcuts ^J How to Add J\ What's New 



Current Folder ■«- p ? x 

I'J « M-File- ► -*~P © *' 



Name j - 



s ..: Archive 
h '..: html 

|_j AcetyleneFJottle.asv 

1^\ AcetyleneFJottle.m 

QProbleml&_Ql.m 

Qprobleml6_02.m 

BProbleml6_03.m 

Bprobleml6_Q4,m 

D Probleml6_D5.asv 

QProbleml6_Q5,m 

□ Probleml6_06.asv 

©Pr€bleml6_Q6,m 

,__, ThermalConductivity.... ■* | 

Acetylene Bottle.m [MATLAB Script) v 

This script computes the 

volume of an acetylene bottle 
with a radius r=0.3 m, a 
hemispherical top 



Command Window ■+■ □ * x 

This is a Classroom License for instructional use only. 
Research and corr.rr.erci a 1 use is prohibited. 

MATE. A3 desktop keyboard shortcuts, such as Ctrl+5, are now customizable. 
In addition, many keyboard shortcuts have changed for improved consistenc 
across the desktop. 

To customize keyboard shortcuts, use Ereferer.ces . From there, you can als 
restore previous default settings by following the steps outlined in Help 



-■L-2'x r.ere :.f you do not ■ 



Vol total = 



to see this me33age again. 



±J 



Workspace 

a <k 1 



BVoLcyl 
HVolJop 
HVolJotal 

h 



Name ■ 



fP Select.,, - 



Value 



0.4241 
0.0565 
0.4807 
1.5000 
0.3000 



Jj 



Command History 



27/09/2011 11:39 
27/09/2011 14:55 



« Start 



Figure 4.4: The MATLAB output in the command window. 



4.1.2 The input Function 

Notice that the script we have created above (Example 4.1) is not interactive and computes the 
total volume only for the variables defined in the m-file. To make this script interactive we will 
make some changes to the existing AcetyleneBottle .m by adding input function and save it as 
AcetyleneBottlelnteractive .m. 

The syntax for input is as follows: 

userResponse = input (' prompt ' ) 

Example 4.2 

Now, let's incorporate the input command in AcetyleneBottlelnteractive .m as 
shown below and the subsequent figure: 
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% This script computes the volume of an acetylene bottle 
% user is prompted to enter 

% a radius r for a hemispherical top 
% a height h for a cylindrical part 
r=input( 'Enter the radius of acetylene bottle in meters '); 

h=input( 'Enter the height of cylindrical part of acetylene bottle in meters '); 
Vol_top=(2*pi*r~3)/3; % Calculating the volume of hemispherical top [m3] 
Vol_cyl=pi*r~2*h; % Calculating the volume of cylindrical bottom [m3] 

Vol_total=Vol_top+Vol_cyl % Calculating the total volume of acetylene bottle [m3] 



agBEHSn&Biaa] 
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SB 
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fflBSD 



-1.0 
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% This script: compares Che volume of ah acetylene bottle 
% user is prompted to enter 

% a radius r for a hemispherical top 

% a height h for a cylindrical part 
r=input ( ' Enter the radius of acetylene bottle in meters ■ ) : 
h=input (■ Enter the height of cylinderical part of acetylene bottle rn meters '); 
Vol_top=(2* , pi*r' - 3t/3; h Calculating Che volume of hemispherical top [n.3] 
Vol_cyl=pi*r~2*h; % Calculating the volume of cylinderical bottom. [ra3] 

Vol total=Vol top+Vol cyl % Calculating Che total volume of acetylene bottle [ 



| script 



In 2 Col 1 



Figure 4.5: Interactive script that computes the volume of acetylene cylinder. 



The command window upon run will be as follows, note that user keys in the radius and 
height values and the same input values result in the same numerical answer as in example 
(Example 4.1) which proves that the computation is correct. 
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Figure 4.6: The same numerical result is obtained through interactive script. 



4.1.3 The disp Function 

As you might have noticed, the output of our script is not displayed in a well-formatted fashion. 
Using disp, we can control how text or arrays are displayed in the command window. For example, 
to display a text string on the screen, type in disp ( 'Hello world! '). This command will return 
our friendly greeting as follows: Hello world! 

disp (variable) can be used to display only the value of a variable. To demonstrate this, issue 
the following command in the command window: 

b = [12 3 4 5] 

We have created a row vector with 5 elements. The following is displayed in the command window: 

> b = [12 3 4 5] 
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Now if we type in disp(b) and press enter, the variable name will not be displayed but its value 
will be printed on the screen: 

> disp(b) 
12 3 4 5 

The following example demonstrates the usage of disp function. 

Example 4.3 

Now, let's open AcetyleneBottlelnteractive .m file and modify it by using the disp 
command. First save the file as AcetyleneBottlelnteractiveDisp.m, so that we don't 
accidentally introduce errors to a working file and also we can easily find this particular 
file that utilizes the disp command in the future. The new file should contain the code 
below: 

% This script computes the volume of an acetylene bottle 

% user is prompted to enter 

% a radius r for a hemispherical top 
% a height h for a cylindrical part 
clc % Clear screen 

disp ('This script computes the volume of an acetylene bottle') 
r=input( 'Enter the radius of acetylene bottle in meters '); 

h=input( 'Enter the height of cylindrical part of acetylene bottle in meters '); 
Vol_top=(2*pi*r~3)/3; % Calculating the volume of hemispherical top [m3] 

Vol_cyl=pi*r~2*h; °/ Calculating the volume of cylindrical bottom [m3] 

Vol_total=Vol_top+Vol_cyl; % Calculating the total volume of acetylene bottle [m3] 

dispC ') % Display blank line 

disp ('The volume of the acetylene bottle is') % Display text 
disp(Vol_total) % Display variable 

Your screen output should look similar to the one below: 

This script computes the volume of an acetylene bottle 
Enter the radius of acetylene bottle in meters .3 
Enter the height of cylindrical part of acetylene bottle in meters 1.5 

The volume of the acetylene bottle is 
0.4807 
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4.1.4 The num2str Function 



The num2str function allows us to convert a number to a text string. Basic syntax is str = 
num2str(A) where variable A is converted to a text and stored in str. Let's see how it works in 
AcetyleneBottlelnteractiveDisp.m. Remember to save the file with a different name before 
editing it, for example, AcetyleneBottlelnteractiveDispl .m. 

Example 4.4 

Add the following line of code to your file: 

str = ['The volume of the acetylene bottle is ', num2str(Vol_total) , ' cubic meters. 

Notice that the three arguments in str are separated with commas. The first argument is 
a simple text that is contained in ' ' . The second argument is where the number to string 
conversion take place. And finally the third argument is also a simple text that completes 
the sentence displayed on the screen. Using semicolon at the end of the line suppresses 
the output. In the next line of our script, we will call str with disp (str) ; . 

AcetyleneBottlelnteractiveDispl.m file should look like this: 

% This script computes the volume of an acetylene bottle 
% user is prompted to enter 

% a radius r for a hemispherical top 

% a height h for a cylindrical part 
clc % Clear screen 

dispCThis script computes the volume of an acetylene bottle:') 
dispC ') % Display blank line 

r=input( 'Enter the radius of acetylene bottle in meters '); 

h=input( 'Enter the height of cylindrical part of acetylene bottle in meters '); 
Vol_top=(2*pi*r~3)/3; % Calculating the volume of hemispherical top [m3] 
Vol_cyl=pi*r~2*h; °/ Calculating the volume of cylindrical bottom [m3] 

Vol_total=Vol_top+Vol_cyl; % Calculating the total volume of acetylene bottle [m3] 
dispC ') % Display blank line 

str = ['The volume of the acetylene bottle is ', num2str(Vol_total) , ' cubic meters.']; 
disp(str) ; 

Running the script should produce the following: 

This script computes the volume of an acetylene bottle: 

Enter the radius of acetylene bottle in meters .3 

Enter the height of cylindrical part of acetylene bottle in meters 1.5 

The volume of the acetylene bottle is 0.48066 cubic meters. 
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4.1.5 The diary Function 

Instead of writing a script from scratch, we sometimes solve problems in the Command Window 
as if we are using a scientific calculator. The steps we perform in this fashion can be used to create 
an m-file. For example, the diary function allows us to record a MATLAB session in a file and 
retrieve it for review. Reviewing the file and by copying relevant parts of it and pasting them in to 
an m-file, a script can be written easily. 

Typing diary at the MATLAB prompt toggles the diary mode on and off. As soon as the diary 
mode is turned on, a file called diary is created in the current directory. If you like to save that 
file with a specific name, say for example probleml6, type diary ( 'probleml6 ' ) . A file named 
problem 16 will be created. The following is the content of a diary file called problem 16. Notice 
that in that session, the user is executing the four files we created earlier. The user's keyboard input 
and the resulting display output is recorded in the file. The session is ended by typing diary which 
is printed in the last line. 

AcetyleneBottle 

Vol.total = 

0.4807 

AcetyleneBottle Interactive 

Enter the radius of acetylene bottle in meters .3 

Enter the height of cylinderical part of acetylene bottle in meters 1.5 

Vol.total = 

0.4807 

AcetyleneBottle Interact iveDisp 

This script computes the volume of an acetylene bottle 

Enter the radius of acetylene bottle in meters .5 

Enter the height of cylinderical part of acetylene bottle in meters 1.6 

The volume of the acetylene bottle is 
1.5184 

AcetyleneBottle Interact iveDispl 

This script computes the volume of an acetylene bottle: 

Enter the radius of acetylene bottle in meters .9 

Enter the height of cylinderical part of acetylene bottle in meters 1.9 
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The volume of the acetylene bottle is 6.3617 cubic meters, 
diary 



4.1.6 Style Guidelines 

Try to apply the following guidelines when writing your scripts: 

• Share your code or programs with others, consider adopting one of Creative Commons 2 or 
GNU General Public License 3 schemes 

• Include your name and contact info in the opening lines 

• Use comments liberally 

• Group your code and use proper indentation 

• Use white space liberally 

• Use descriptive names for your variables 

• Use descriptive names for your m-files 



4.1.7 Summary of Key Points 

1. A script is a file containing a sequence of MATLAB statements. Script files have a filename 
extension of .m. 

2. Functions such as input, disp and num2str can be used to make scripts interactive, 

3. diary function is useful to record a MATLAB command window session from which an 
m-file can be easily created, 

4. Various style guidelines covered here help improve our code. 



4.2 Problem Set 4 

Exercise 4.2.1 (Solution on p. 100.) 

Write a script that will ask for pressure value in psi and display the equivalent pressure in 
kPa with a statement, such as "The converted pressure is: ..." 

Exercise 4.2.2 (Solution on p. 100.) 

Write a script that generates a table of conversions from Fahrenheit to Celsius tempera- 
tures for a range and increment entered by the user, such as 

Enter the beginning temperature in F: 
Enter the ending temperature in F: 



2 http://creativecommons.org/ 
3 http://www.gnu.org/licenses/gpl-3.0.html 

4n 



This content is available online at <http://cnx.org/content/m41536/!. 2/>. 
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Enter the increment value: 

Test your script with 20 the beginning Fahrenheit value, 200 the ending Fahrenheit 
value and 20 the increment. 

Exercise 4.2.3 (Solution on p. 101.) 

Pascal's Law states that pressure is transmitted undiminished in all directions throughout 
a fluid at rest. (See the illustration below). An initial force of 150 N is transmitted from 
a piston of 25 mm A 2 to a piston of 100 mm A 2. This force is progressively increased up 
to 200 N. Write a script that computes the corresponding load carried by the larger piston 
and tabulate your results. 



A 1 



P1 



P2 



A 2 



Figure 4.7: A simple hydraulic system. 



Exercise 4.2.4 (Solution on p. 102.) 

Modify your script in previous problem (Exercise 4.2.3) so that the user provides the 
following input: 



Enter the initial force in N: 

Enter the final force in N: 

Enter the increment value: 

Enter the area of small piston in mm A 2: 

Enter the area of big piston in mm A 2: 
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Test your script with 150, 200, 10, 25 and 100 with respect to each input variable. 

Exercise 4.2.5 (Solution on p. 103.) 

Write a script to solve the Stress-Strain problem in the Problem Set (Problem 2.2.9) 

Exercise 4.2.6 (Solution on p. 104.) 

Modify the script, you wrote above (Exercise 4.2.5)and plot an annotated Stress-Strain 
graph. 
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Solutions to Exercises in Chapter 4 

Solution to Exercise 4.2.1 (p. 97) 



% This script converts pressures from psi to kPa 

% User is prompted to enter pressure in psi 

clc % Clear screen 

dispCThis script converts pressures from psi to kPa:') 

dispC ') % Display blank line 

psi=input('What is the pressure value in psi? '); 

kPa=psi*6. 894757; % Calculating pressure in kPa 

dispC ') % Display blank line 

str = ['The converted pressure is: ', num2str(kPa) , ' kPa. '] ; 

disp(str) ; 

The script output is as follows: 



This script converts pressures from psi to kPa: 

What is the pressure value in psi? 150 

The converted pressure is: 1034.2135 kPa. 
Solution to Exercise 4.2.2 (p. 97) 

'/, This script generates a table of conversions 
% From Fahrenheit to Celsius temperatures 
clc % Clear screen 

dispCThis script generates a table of conversions from Fahrenheit to Celsius') 
dispC ') % Display blank line 

lowerF=input( 'Enter the beginning temperature in F: '); 
upperF=input( 'Enter the ending temperature in F: '); 
increment=input( 'Enter the increment value: '); 

Fahrenheit= [lowerF: increment :upperF] ; °/„ Creating a row vector with F values 
Celsius=5/9*(Fahrenheit-32) ; °/ Converting from F to C 
dispC ') % Display blank line 

str = ['Fahrenheit Celsius '] ;°/ Displaying table header 
disp(str) ; 

% Tabulating results in two columns, ' is being used to transpose row to column 
disp( [Fahrenheit' Celsius']) 
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The script output is as follows: 

This script generates a table of conversions from Fahrenheit to Celsius 

Enter the beginning temperature in F: 20 
Enter the ending temperature in F: 200 
Enter the increment value : 20 

Fahrenheit Celsius 

20.0000 -6.6667 

40.0000 4.4444 

60.0000 15.5556 

80.0000 26.6667 

100.0000 37.7778 

120.0000 48.8889 

140.0000 60.0000 

160.0000 71.1111 

180.0000 82.2222 

200.0000 93.3333 

Solution to Exercise 4.2.3 (p. 98) 



% This script computes the load carried by the larger piston in a hydraulic system 

clc % Clear screen 

dispCThis script computes the load carried by the larger piston in a hydraulic system') 

dispC ') % Display blank line 

initialF=150; 

finalF=200; 

increment=10; 

areal=25; 

area2=100; 

Fl=[initialF: increment :finalF] ; °/„ Creating a row vector with Fl values 

F2=Fl*area2/areal; °/„ Calculating F2 values 

dispC ') % Display blank line 

str = [' Fl F2 '];% Displaying table header 

disp(str) ; 

disp([Fl' F2']) °/ Tabulating results in two columns, ' is being used to transpose row to 

The script output is as follows: 



102 CHAPTER 4. INTRODUCTORY PROGRAMMING 

This script computes the load carried by the larger piston in a hydraulic system 



Fl F2 

150 600 

160 640 

170 680 

180 720 

190 760 

200 800 



Solution to Exercise 4.2.4 (p. 98) 



% This script computes the load carried by the larger piston in a hydraulic system 

clc % Clear screen 

dispCThis script computes the load carried by the larger piston in a hydraulic system') 

dispC ') % Display blank line 

initialF=input( 'Enter the initial force in N: '); 

f inalF= input ( 'Enter the final force in N: '); 

increment=input( 'Enter the increment value: '); 

areal=input(' Enter the area of small piston in mm~2: '); 

area2=input( 'Enter the area of big piston in mm~2: '); 

Fl=[initialF: increment :finalF] ; °/„ Creating a row vector with Fl values 

F2=Fl*area2/areal; °/„ Calculating F2 values 

dispC ') % Display blank line 

str = [' Fl F2 '];% Displaying table header 

disp(str) ; 

disp([Fl' F2']) °/ Tabulating results in two columns, ' is being used to transpose row to 

The script output is as follows: 

This script computes the load carried by the larger piston in a hydraulic system 

Enter the initial force in N: 150 

Enter the final force in N: 200 

Enter the increment value: 10 

Enter the area of small piston in mm~2: 25 

Enter the area of big piston in mm~2: 100 

Fl F2 
150 600 
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160 640 

170 680 

180 720 

190 760 

200 800 

Solution to Exercise 4.2.5 (p. 99) 

The m-file contains the following: 

% This script uses readings from a Tensile test and 
% Computes Strain and Stress values 
clc % Clear screen 

dispCThis script uses readings from a Tensile test and') 
disp( 'Computes Strain and Stress values') 
dispC ') % Display a blank line 

Specimen_dia=12.7; % Specimen diameter in mm 
*/. Load in kM 
Load_kN=[0;4.89;9.779;14.67;19.56;24.45; . . . 

27. 62; 29. 39; 32. 68 ;33. 95; 34. 58; 35. 22 ;.. . 

35. 72; 40. 54; 48. 39; 59. 03; 65. 87; 69. 42 ;.. . 

69.67;68.15;60.81]; 
% Gage length in mm 
Length_mm= [50.8; 50 . 8102 ; 50 . 8203 ; 50 . 8305 ; . . . 

50.8406;50.8508;50.8610;50.8711; . . . 

50.9016;50.9270;50.9524;50.9778; . . . 

51.0032;51.816;53.340;55.880;58.420; . . . 

60.96;61.468;63.5;66.04] ; 

% Calculate x-sectional area im m2 

Cross_sectional_Area=pi/4* ( (Specimen_dia/1000) "2) ; 

% Calculate change in length, initial lenght is 50.8 mm 

Delta_L=Length_mm-50 . 8 ; 

% Calculate Stress in MPa 

Sigma=(Load_kN./Cross_sectional_Area)*10~(-3) ; 

% Calculate Strain in mm/mm 

Epsilon=Delta_L . /50 . 8 ; 

str = ['Specimen diameter is ', num2str(Specimen_dia) , ' mm.']; 

disp(str) ; 

Results=[Load_kN Length_mm Delta_L Sigma Epsilon] ; 

% Tabulated results 

dispC Load Length Delta L Stress Strain') 

disp(Results) 

After executed, the command window output is: 
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This script uses readings from a Tensile test and 
Computes Strain and Stress values 



Specimen diameter is 12.7 mm. 





Load 


Length 


Delta L 


Stress 


Strain 







50 


8000 

















4 


8900 


50 


8102 





0102 


38 


6022 





0002 


9 


7790 


50 


8203 





0203 


77 


1964 





0004 


14 


6700 


50 


8305 





0305 


115 


8065 





0006 


19 


5600 


50 


8406 





0406 


154 


4086 





0008 


24 


4500 


50 


8508 





0508 


193 


0108 





0010 


27 


6200 


50 


8610 





0610 


218 


0351 





0012 


29 


3900 


50 


8711 





0711 


232 


0076 





0014 


32 


6800 


50 


9016 





1016 


257 


9792 





0020 


33 


9500 


50 


9270 





1270 


268 


0047 





0025 


34 


5800 


50 


9524 





1524 


272 


9780 





0030 


35 


2200 


50 


9778 





1778 


278 


0302 





0035 


35 


7200 


51 


0032 





2032 


281 


9773 





0040 


40 


5400 


51 


8160 


1 


0160 


320 


0269 





0200 


48 


3900 


53 


3400 


2 


5400 


381 


9955 





0500 


59 


0300 


55 


8800 


5 


0800 


465 


9888 





1000 


65 


8700 


58 


4200 


7 


6200 


519 


9844 





1500 


69 


4200 


60 


9600 


10 


1600 


548 


0085 





2000 


69 


6700 


61 


4680 


10 


6680 


549 


9820 





2100 


68 


1500 


63 


5000 


12 


7000 


537 


9830 





2500 


60 


8100 


66 


0400 


15 


2400 


480 


0403 





3000 



Solution to Exercise 4.2.6 (p. 99) 

Edited script contains the plot commands: 

% This script uses readings from a Tensile test and 

% Computes Strain and Stress values 

clc % Clear screen 

dispCThis script uses readings from a Tensile test and') 

disp( 'Computes Strain and Stress values') 

dispC ') % Display a blank line 

Specimen_dia=12.7; % Specimen diameter in mm 

*/. Load in kM 

Load_kN= [0 ; 4 . 89 ; 9 . 779 ; 14 . 67 ; 19 . 56 ; 24 . 45 : 

27. 62 ; 29. 39 ; 32. 68 ;33. 95; 34. 58; 35. 22: 

35. 72; 40. 54; 48. 39; 59. 03; 65. 87; 69. 42: 
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69.67;68.15;60.81]; 
% Gage length in mm 
Length_mm= [50.8; 50 . 8102 ; 50 . 8203 ; 50 . 8305 ; . . . 

50.8406;50.8508;50.8610;50.8711; . . . 

50.9016;50.9270;50.9524;50.9778; . . . 

51.0032;51.816;53.340;55.880;58.420; . . . 

60.96;61.468;63.5;66.04] ; 

% Calculate x-sectional area im m2 

Cross_sectional_Area=pi/4* ( (Specimen_dia/1000) "2) ; 

% Calculate change in length, initial lenght is 50.8 mm 

Delta_L=Length_mm-50 . 8 ; 

% Calculate Stress in MPa 

Sigma=(Load_kN./Cross_sectional_Area)*10~(-3) ; 

% Calculate Strain in mm/mm 

Epsilon=Delta_L . /50 . 8 ; 

str = ['Specimen diameter is ', num2str(Specimen_dia) , ' mm.']; 

disp(str) ; 

Results=[Load_kN Length_mm Delta_L Sigma Epsilon] ; 

% Tabulated results 

dispC Load Length Delta L Stress Strain') 

disp(Results) 

% Plot Stress versus Strain 

plot (Epsilon, Sigma) 

title ('Stress versus Strain Curve') 

xlabel ( ' Strain [mm/mm] ' ) 

ylabelC Stress [mPa]') 

grid 

In addition to Command Window output, the following plot is generated: 
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Stress versus Strain Curve 



Cl- 



io 
m 

W 



600 



500 



400 



300 



200 



100 -■ 



0.05 0.1 



0.15 0.2 

Strain [mm/mm] 



- i r i i r i ■ 



0.25 0.3 0.35 



Chapter 5 
Interpolation 



5.1 Interpolation 




Linear interpolation is one of the most common techniques for estimating values between two 
given data points. For example, when using steam tables, we often have to carry out interpolations. 
With this technique, we assume that the function between the two points is linear. MATLAB has a 



This content is available online at <http://cnx.org/content/m4 1455/1. 2/>. 
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built-in interpolation function. 

5.1.1 The interpl Function 

Give an x-y table, y_new = interpl (x,y,x_new) interpolates to find y_new. Consider the fol- 
lowing examples: 

Example 5.1 

To demonstrate how the interpl function works, let us create an x-y table with the fol- 
lowing commands; 

x = 0:5; 

y = [0,20,60,68,77,110] ; 

To tabulate the output, we can use 

Cx',y'] 

The result is 

ans = 









1 


20 


2 


60 


3 


68 


4 


77 


5 


110 



Suppose we want to find the corresponding value for 1.5 or interpolate for 1.5. Using 
y_new = interpl (x,Y,x_new) syntax, let us type in: 

y_new=interpl (x , y , 1 . 5) 
y_new = 

40 

Example 5.2 

The table we created above has only 6 elements in it and suppose we need a more detailed 
table. In order to do that, instead of a single new x value, we can define an array of new x 
values, the interpl function returns an array of new y values: 
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new_x = 0:0.2:5; 

new_y = interpl(x,y,new_x) ; 

Let's display this table 



[new_x' ,new_y'] 
The result is 

ans = 













0. 


.2000 


4. 


.0000 


0. 


.4000 


8. 


.0000 


0. 


.6000 


12. 


.0000 


0. 


.8000 


16. 


.0000 


1. 


.0000 


20 


.0000 


1. 


.2000 


28 


.0000 


1. 


.4000 


36 


.0000 


1. 


.6000 


44 


.0000 


1. 


.8000 


52 


.0000 


2 


.0000 


60 


.0000 


2 


.2000 


61 


.6000 


2 


.4000 


63 


.2000 


2 


.6000 


64 


.8000 


2 


.8000 


66. 


.4000 


3 


.0000 


68. 


.0000 


3 


.2000 


69 


.8000 


3 


.4000 


71. 


.6000 


3 


.6000 


73 


.4000 


3 


.8000 


75 


.2000 


4. 


.0000 


77 


.0000 


4. 


.2000 


83. 


.6000 


4. 


.4000 


90. 


.2000 


4. 


.6000 


96. 


.8000 


4. 


.8000 


103. 


.4000 


5. 


.0000 


110. 


.0000 



Example 5.3 

Using the table below, find the internal energy of steam at 215 °C and the temperature if 
the internal energy is 2600 kJ/kg (use linear interpolation). 
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Temperature [C] 


Internal Energy [kj/kg] 


100 


2506.7 


150 


2582.8 


200 


2658.1 


250 


2733.7 


300 


2810.4 


400 


2967.9 


500 


3131.6 



Table 5.1: An extract from Steam Tables 
First let us enter the temperature and energy values 

temperature = [100, 150, 200, 250, 300, 400, 500]; 
energy = [2506.7, 2582.8, 2658.1, 2733.7, 2810.4, 2967.9, 3131.6]; 
[temperature ' , energy ' ] 

returns 



ans 



1.0e+003 * 






0.1000 


2 


5067 


0.1500 


2 


5828 


0.2000 


2 


6581 


0.2500 


2 


7337 


0.3000 


2 


8104 


0.4000 


2 


9679 


0.5000 


3 


1316 



issue the following for the first question: 

new_energy = interpl (temperature, energy, 215) 
returns 

new_energy = 



2.6808e+003 



Ill 



Now, type in the following for the second question: 

new_temperature = interpl (energy, temperature, 2600) 
returns 

new_temperature = 
161.4210 



5.1.2 Summary of Key Points 

1. Linear interpolation is a technique for estimating values between two given data points, 

2. Problems involving steam tables often require interpolated data, 

3. MATLAB has a built-in interpolation function. 



5.2 Problem Set 2 

Exercise 5.2.1 (Solution on p. 114.) 

Determine the saturation temperature, specific liquid enthalpy, specific enthalpy of evap- 
oration and specific enthalpy of dry steam at a pressure of 2.04 MPa. 



Pressure [MN/m 2 ] 


Saturation Temperature [C] 


h f [kj/kg] 


h fg [kj/kg] 


h g [kj/kg] 


2.1 


214.9 


920.0 


1878.2 


2798.2 


2.0 


212.4 


908.6 


1888.6 


2797.2 



Table 5.2: An extract from steam tables 



Exercise 5.2.2 (Solution on p. 114.) 

The following table gives data for the specific heat as it changes with temperature for a 
perfect gas. 3 



2 This content is available online at <http://cnx.org/content/m41624/Ll/>. 

3 Thermodynamics and Heat Power by Kurt C. Rolle, Pearson Prentice Hall. ©2005, (p. 19) 
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Temperature [F] 


Specific Heat [BTU/lbmF] 


25 


0.118 


50 


0.120 


75 


0.123 


100 


0.125 


125 


0.128 


150 


0.131 



Table 5.3: Change of specific heat with temperature 

Using interpl function calculate the specific heat for 30 F, 70 F and 145 F. 

Exercise 5.2.3 (Solution on p. 115.) 

For the problem above (Exercise 5.2.2), create a more detailed table in which temperature 
varies between 25 and 150 with 5 F increments and corresponding specific heat values. 

Exercise 5.2.4 (Solution on p. 116.) 

During a 12-hour shift a fuel tank has varying levels due to consumption and transfer 
pump automatically cutting in and out to maintain a safe fuel level. The following table of 
fuel tank level versus time is missing readings for 5 and 9 AM. Using linear interpolation, 
estimate the fuel level at those times. 



Time [hours, AM] 


Tank level [m] 


1:00 


1.5 


2:00 


1.7 


3:00 


2.3 


4:00 


2.9 


5:00 


? 


6:00 


2.6 


7:00 


2.5 


8:00 


2.3 


9:00 


? 


10:00 


2.0 


11:00 


1.8 


12:00 


1.3 
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Table 5.4: Fuel tank level versus time 
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Solutions to Exercises in Chapter 5 

Solution to Exercise 5.2.1 (p. Ill) 

MATLAB solution is as follows; 

> pressure=[2.1 2.0] ; 

> sat_temp=[214.9 212.4]; 

> h_f=[920 908.6]; 

> h_fg=[1878.2 1888.6]; 

> h_g=[2798.2 2797.2] ; 

3> sat_temp_new=interpl (pressure , sat_temp , 2 . 04) 
sat_temp_new = 

213.4000 
» h_f _new=interpl (pressure , h_f , 2 . 04) 
h_f_new = 

913.1600 
3> h_f g_new=interpl (pressure ,h_fg, 2 . 04) 
h_fg_new = 

1.8844e+003 
» h_g_new=interpl (pressure , h_g , 2 . 04) 
h_g_new = 

2.7976e+003 

Solution to Exercise 5.2.2 (p. Ill) 

MATLAB solution is as follows: 

> temperature= [25 ; 50 ; 75 ; 100 ; 125 ; 150] 
temperature = 

25 



115 



50 

75 
100 
125 
150 

> specif ic_heat=[. 118; .120; .123; .125; .128; .131] 

specific_heat = 

0.1180 
0.1200 
0.1230 
0.1250 
0.1280 
0.1310 

» specif ic_heat At30=interpl (temperature , specif ic_heat , 30) 

specific_heatAt30 = 

0.1184 
» specif ic_heat At70=interpl (temperature , specif ic_heat , 70) 
specific_heatAt70 = 

0.1224 
» specif ic_heatAtl45=interpl (temperature , specif ic_heat , 145) 
specific_heatAtl45 = 

0.1304 



Solution to Exercise 5.2.3 (p. 112) 

MATLAB solution is as follows: 

3> new_temperature=25:5: 150; 
^> new_specif ic_heat=interpl (temperature , specif ic_heat ,new_temperature) ; 
» [new_temperature' ,new_specif ic_heat '] 
ans = 
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25.0000 0.1180 

30.0000 0.1184 

35.0000 0.1188 

40.0000 0.1192 

45.0000 0.1196 

50.0000 0.1200 

55.0000 0.1206 

60.0000 0.1212 

65.0000 0.1218 

70.0000 0.1224 

75.0000 0.1230 

80.0000 0.1234 

85.0000 0.1238 

90.0000 0.1242 

95.0000 0.1246 

100.0000 0.1250 

105.0000 0.1256 

110.0000 0.1262 

115.0000 0.1268 

120.0000 0.1274 

125.0000 0.1280 

130.0000 0.1286 

135.0000 0.1292 

140.0000 0.1298 

145.0000 0.1304 

150.0000 0.1310 

Solution to Exercise 5.2.4 (p. 112) 

> time=[l 2 3 4 6 7 8 10 11 12] ; 
> tank_level=[1.5 1.7 2.3 2.9 2.6 2.5 2.3 2.0 1.8 1.3]; 

3> tank_level_at_5=interpl (time , tank_level , 5) 

tank_level_at_5 = 

2.7500 

3> tank_level_at_9=interpl (time , tank_level , 9) 
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tank_level_at_9 = 
2.1500 
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Chapter 6 
Numerical Integration 



6.1 Computing the Area Under a Curve 1 




This chapter essentially deals with the problem of computing the area under a curve. First, we 
will employ a basic approach and form trapezoids under a curve. From these trapezoids, we can 
calculate the total area under a given curve. This method can be tedious and is prone to errors, so in 



This content is available online at <http://cnx.org/content/m4 1454/1. 3/>. 
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the second half of the chapter, we will utilize a built-in MATLAB function to carry out numerical 
integration. 

6.1.1 A Basic Approach 

There are various methods to calculating the area under a curve, for example, Rectangle Method 2 
, Trapezoidal Rule 3 and Simpson's Rule 4 . The following procedure is a simplified method. 

Consider the curve below: 



7o 



7i 


y-i 


y-i 


y* 


y$ 



AX AX AX AX AX 



Figure 6.1: Numerical integration 



Each segment under the curve can be calculated as follows: 

- (y + yi) Ax + - (yi +y 2 ) Ax + - (y 2 +y3) Ax (6.1) 

Therefore, if we take the sum of the area of each trapezoid, given the limits, we calculate the total 
area under a curve. Consider the following example. 

Example 6.1 

Given the following data, plot an x-y graph and determine the area under a curve between 

x=3 and x=30 



2 http://en. wikipedia.org/wiki/Rectangle_method 
3 http://en. wikipedia.org/wiki/Trapezoidal_rule 
4 http://en. wikipedia.org/wiki/Simpson%27s_rule 
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Index 


x[m] 


y[N] 


1 


3 


27.00 


2 


10 


14.50 


3 


15 


9.40 


4 


20 


6.70 


5 


25 


5.30 


6 


30 


4.50 



Table 6.1: Data Set 

First, let us enter the data set. For x, issue the following command 

x=[3, 10, 15,20,25,30] ; . And for y, y=[27, 14.5,9.4,6.7,5.3,4.5] ;. If yu type 
in [x ' , y ' ] , you will see the following tabulated result. Here we transpose row vectors 
with ' and displaying them as columns: 

ans = 



3 


0000 


27 


0000 


10 


0000 


14 


5000 


15 


0000 


9 


4000 


20 


0000 


6 


7000 


25 


0000 


5 


3000 


30 


0000 


4 


5000 



Compare the data set above with the given information in the question (Table 6.1). 
To plot the data type the following: 

plot(x,y) ,title( 'Distance-Force Graph') , xlabelC Distance [m] ') , ylabelC Force [N] ') ,gri 
The following figure is generated: 
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30 



Distance-Force Graph 



20 



15 



ID 






ID 15 

Distancefm] 



20 



30 



Figure 6.2: Distance-Force Graph 



To compute dx for consecutive x values, we will use the index for each x value, see the 
given data in the question (Table 6.1).: 

dx=[x(2)-x(l) J x(3)-x(2) J x(4)-x(3) J x(5)-x(4) J x(6)-x(5)]; 

dy is computed by the following command: 

dy=[0.5*(y(2)+y(l)) J 0.5*(y(3)+y(2)) 5 0.5*(y(4)+y(3)) J 0.5*(y(5)+y(4)) J 0.5*(y(6)+y(5))] 

dx and dy can be displayed with the following command: [dx ' , dy ' ] . The result will look 
like this: 

[dx'.dy'] 



ans 










7.0000 


20 


7500 




5.0000 


11 


9500 




5.0000 


8 


0500 




5.0000 


6 


0000 




5.0000 


4 


9000 
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Our results so far are shown below 



x[m] 


y[N] 


dx [m] 


dy[N] 


3 


27.00 






10 


14.50 


7.00 


20.75 


15 


9.40 


5.00 


11.95 


20 


6.70 


5.00 


8.05 


25 


5.30 


5.00 


6.00 


30 


4.50 


5.00 


4.90 



Table 6.2: x, y and corresponding differential elements 

If we multiply dx by dy, we find da for each element under the curve. The differential 
area da=dx*dy, can be computed using the 'term by term multiplication' technique in 
MATLAB as follows: 

da=dx . *dy 



da 



145.2500 59.7500 40.2500 30.0000 24.5000 

Each value above represents an element under the curve or the area of trapezoid. By taking 
the sum of array elements, we find the total area under the curve. 

sum (da) 



ans 



299.7500 
The following (Table 6.3) illustrates all the steps and results of our MATLAB computation. 



x[m] 


y[N] 


dx [m] 


dy[N] 


dA [Nm] 


3 


27.00 








10 


14.50 


7.00 


20.75 


145.25 


15 


9.40 


5.00 


11.95 


59.75 


20 


6.70 


5.00 


8.05 


40.25 


25 


5.30 


5.00 


6.00 


30.00 


30 


4.50 


5.00 


4.90 


24.50 










299.75 
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Table 6.3: Computation of the approximate area under a curve 

6.1.2 The Trapezoidal Rule 

Sometimes it is rather convenient to use a numerical approach to solve a definite integral. The 
trapezoid rule allows us to approximate a definite integral using trapezoids. 

6.1.2.1 The trapz Command 

Z = trapz(Y) computes an approximation of the integral of Y using the trapezoidal method. 

Now, let us see a typical problem. 
Example 6.2 

Given Area = f 2 x 2 dx, an analytical solution would produce 39. Use trapz command and 
solve it 

1 . Initialize variable x as a row vector, from 2 with increments of 0. 1 to 5: x=2 : .1:5; 

2. Declare variable y as y=x~2;. Note the following error prompt: ??? Error 
using ==> mpower Inputs must be a scalar and a square matrix. 
This is because x is a vector quantity and MATLAB is expecting a scalar input for 
y. Because of that, we need to compute y as a vector and to do that we will use the 
dot operator as follows: y=x . "2 ; . This tells MATLAB to create vector y by taking 
each x value and raising its power to 2. 

3. Now we can issue the following command to calculate the first area, the output will 
be as follows: 

areal=trapz (x , y) 
areal = 

39.0050 

Notice that this numerical value is slightly off. So let us increase the number of increments 
and calculate the area again: 

x=2: .01:5; 
y=x.~2; 
area2=trapz (x , y) 

area2 = 

39.0001 
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Yet another increase in the number of increments: 

x=2: .001:5; 
y=x.~2; 
area3=trapz (x , y) 

area3 = 

39.0000 

Example 6.3 

Determine the value of the following integral: 

Jq sin (x) dx 

1. Initialize variable x as a row vector, from with increments of pi/100 to pi: 
x=0:pi/100:pi; 

2. Declare variable y as y=sin(x) ; 

3. Issue the following command to calculate the first area, the output will be as follows: 

areal=trapz (x , y) 

areal = 

1.9998 

let us increase the increments as above: 

x=0:pi/1000:pi; 
y=sin(x) ; 
area2=trapz (x , y) 

area2 = 

2.0000 



6.1.3 Summary of Key Points 

1 . In its simplest form, numerical integration involves calculating the areas of segments that 
make up the area under a curve, 

2. MATLAB has built-in functions to perform numerical integration, 

3. Z = trapz(Y) computes an approximation of the integral of Y using the trapezoidal 
method. 
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6.2 Problem Set 5 

Exercise 6.2.1 (Solution on p. 128.) 

Let the function y defined by y = cos (x). Plot this function over the interval [-pi,pi]. Use 
numerical integration techniques to estimate the integral of y over [0, pi] and over [-pi, pi]. 

Exercise 6.2.2 (Solution on p. 128.) 

Let the function y defined by y — 0.04x 2 — 2.1 3x + 32. 58. Plot this function over the 
interval [3,30]. Use numerical integration techniques to estimate the integral of y over 
[3,30]. 

Exercise 6.2.3 (Solution on p. 129.) 

A 2000-liter tank is full of lube oil. It is known that if lube oil is drained from the tank, 
the mass flow rate will decrease from the maximum when the tank level is at the highest. 
The following data were collected when the tank was drained. 



Time [min] 


Mass Flow [kg/min] 





50.00 


5 


48.25 


10 


46.00 


15 


42.50 


20 


37.50 


25 


30.50 


30 


19.00 


35 


9.00 



Table 6.4: Data 



Write a script to estimate the amount of oil drained in 35 minutes. 

Exercise 6.2.4 (Solution on p. 130.) 

A gas is expanded in an engine cylinder, following the law PV 1 3 =c. The initial pressure 
is 2550 kPa and the final pressure is 210 kPa. If the volume at the end of expansion is 0.75 
m 3 , compute the work done by the gas. 6 

Exercise 6.2.5 (Solution on p. 130.) 

A force F acting on a body at a distance s from a fixed point is given by F — 3s + \ . Write 
a script to compute the work done when the body moves from the position where s=l to 
that where s=10. 7 



5 This content is available online at <http://cnx.org/content/m41541/L7/>. 

6 Applied Heat for Engineers by W. Embleton and L Jackson, Thomas Reed Publications. ©1999, (p. 80) 

1 0. N. Mathematics: 2 by J. Dobinson, Penguin Library of Technology. ©1969, (p. 213) 
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Exercise 6.2.6 (Solution on p. 131.) 

The pressure p and volume v of a given mass of gas are connected by the relation 

[P + "T ) ( v ~~ b) — k where a, b and k are constants. Express p in terms of v, and write a 
script to compute the work done by the gas in expanding from an initial volume to a final 
volume. 8 

Test your solution with the following input: 

a: 0.01 

b: 0.001 

The initial pressure [kPa]: 100 

The initial volume [m3]: 1 

The final volume [m3]: 2 



O. N. Mathematics: 2 by J. Dobinson, Penguin Library of Technology. ©1969, (p. 212) 
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Solutions to Exercises in Chapter 6 

Solution to Exercise 6.2.1 (p. 126) 

1. Plotting: 

x=-pi :pi/100:pi; 

y=cos(x) ; 

plot (x , y) , title ( ' Graph of y=cos (x) ' ) , xlabel ( ' x ' ) , ylabel ( ' y ' ) , grid 

2. Area calculation 1 : 

> x=0:pi/100:pi; 

> y=cos(x) ; 

^> areal=trapz(x,y) 

areal = 

1.1796e-016 

3. Area calculation 2: 

> x=-pi:pi/100:pi; 

> y=cos(x) ; 

3> area2=trapz(x,y) 

area2 = 

1.5266e-016 
Solution to Exercise 6.2.2 (p. 126) 

1. Plotting: 

> x=3: .1:30; 

> y=0.04*(x.~2)-2.13.*x+32.58; 
3> plot(x,y), titleCGraph of . . . 
y=.04*(x~2)-2.13*x+32.58'),xlabel('x'),ylabel('y'),grid 

2. Area calculation: 

^> area=trapz(x,y) 
area = 
290.3868 
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Solution to Exercise 6.2.3 (p. 126) 

clc 
t=linspace (0,35,8) % Data entry for time [min] 

m=[50 48.25 46 42.5 37.5 30.5 19 9] °/„ Data entry for mass flow [kg/min] 
% Calculate time intervals 
dt=[t(2)-t(l) ,t(3)-t(2) ,t(4)-t(3) , . . . 
t(5)-t(4),t(6)-t(5),t(7)-t(6),t(8)-t(7)] 
% Calculate mass out 

dm=[0.5*(m(2)+m(l)),0.5*(m(3)+m(2)),0.5*(m(4)+m(3)),0.5*(m(5)+. . . 
m(4)) ,0.5*(m(6)+m(5)) ,0.5*(m(7)+m(6)) ,0.5*(m(8)+m(7))] 
% Calculate differential areas 
da=dt . *dm; 

°/. Tabulate time and mass flow 
Ct',m'] 

% Tabulate time intervals, mass out and differential areas 
[dt',dm',da'] 

% Calculate the amount of oil drained [kg] in 35 minutes 
Oil_Drained=sum (da) 

The output is: 



ans 








50. 


.0000 


5. 


.0000 


48 


.2500 


10 


.0000 


46 


.0000 


15 


.0000 


42 


.5000 


20 


.0000 


37 


.5000 


25 


.0000 


30 


.5000 


30 


.0000 


19 


.0000 


35. 


.0000 


9. 


.0000 



ans 



5.0000 49.1250 245.6250 

5.0000 47.1250 235.6250 

5.0000 44.2500 221.2500 

5.0000 40.0000 200.0000 

5.0000 34.0000 170.0000 

5.0000 24.7500 123.7500 



130 



CHAPTER 6. NUMERICAL INTEGRATION 



5.0000 14.0000 70.0000 



Oil_Drained = 

1.2663e+003 
Solution to Exercise 6.2.4 (p. 126) 

clc 

disp('A gas is expanded in an engine cylinder, following the law PV~1.3=c') 

dispCThe initial pressure is 2550 kPa and the final pressure is 210 kPa.') 

dispClf the volume at the end of expansion is 0.75 m3,') 

dispCCompute the work done by the gas.') 

dispC ') % Display blank line 

n=1.3; 

P_i=2550; % Initial pressure 

P_f=210; °/„ Final pressure 

V_f=.75; % Final volume 

V_i=(P_f*(V_f~n)/P_i)~(l/n); °/„ Initial volume 

c=P_f*V_f~n; 

v=V_i:.001:V_f; 

p=c./(v.~n) ; 

WorkDone=trapz ( v , p) 



% Creating a row vector for volume, v 
% Computing pressure for volume 
% Integrating p*dv 



The output is: 



A gas is expanded in an engine cylinder, following the law PV~1.3=c 

The initial pressure is 2550 kPa and the final pressure is 210 kPa. 

If the volume at the end of expansion is 0.75 m3, 
Compute the work done by the gas . 



WorkDone = 



409.0666 
Solution to Exercise 6.2.5 (p. 126) 



clc 
disp('A force F acting on a body at a distance s from a fixed point is given by') 
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disp('F=3*s+(l/(s~2)) where s is the distance in meters') 
disp( 'Compute the total work done in moving') 
dispCFrom the position where s=l to that where s=10.') 
dispC ') % Display blank line 

s=l:. 001:10; % Creating a row vector for distance, s 
F=3.*s+(l./(s.~2)) ; °/„ Computing Force for s 
WorkDone=trapz(s,F) °/ Integrating F*ds over 1 to 10 meters. 

The output is: 

A force F acting on a body at a distance s from a fixed point is given by 
F=3*s+(l/(s~2)) where s is the distance in meters 
Compute the total work done in moving 
From the position where s=l to that where s=10. 



WorkDone = 

149.4000 
Solution to Exercise 6.2.6 (p. 126) 

clc % Clear screen 

dispCThis script computes the work done by') 
dispCThe gas in expanding from volume vl to v2') 
dispC ') % Display blank line 

a=input( 'Enter the constant a: '); 
b=input( 'Enter the constant b: '); 
p_i=input( 'Enter the initial pressure [kPa] : '); 
v_i=input( 'Enter the initial volume [m3] : '); 
v_f=input( 'Enter the final volume [m3] : '); 
k=(p_i+(a/(v_i~2))*(v_i-b)) ; % Calculating constant k 
v=v_i : .001 :v_f ; °/„ Creating a row vector for volume 

p=(k./(v-b))-(a./(v.~2)) ; % Computing pressure for volume 
WorkDone=trapz(v,p) ; °/„ Integrating p*dv 
dispC ') % Display blank line 

str = ['The work done by the gas in expanding from ', num2str(v_i) , 

' m3 to ' num2str(v_f) , ' m3 is ', num2str (WorkDone) , ' kW. '] ; 
disp(str) ; 

The output is: 

This script computes the work done by 
The gas in expanding from volume vl to v2 
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Enter the constant a: .01 

Enter the constant b: .001 

Enter the initial pressure [kPa] : 100 

Enter the initial volume [m3] : 1 

Enter the final volume [m3] : 2 

The work done by the gas in expanding from 1 m3 to 2 m3 is 69.3667 kW. 



Chapter 7 
Regression Analysis 



7.1 Linear Regression 




'This content is available online at <http://cnx.org/content/m4 1448/1 .2/>. 
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7.1.1 What is Regression Analysis? 

Suppose we calculate some variable of interest, y, as a function of some other variable x. We 
call y the dependent variable and x the independent variable. For example, consider the data set 
below, taken from a simple experiment involving a vehicle, its velocity versus time is tabulated. In 
this case, velocity is a function of time, thus velocity is the dependent variable and the time is the 
independent variable. 



Time [s] 


Velocity [m/s] 





20 


10 


39 


20 


67 


30 


89 


40 


111 


50 


134 


60 


164 


70 


180 


80 


200 



Table 7.1: Vehicle velocity versus time. 

In its simplest form regression analysis involves fitting the best straight line relationship to explain 
how the variation in a dependent variable, y, depends on the variation in an independent variable, x. 
In our example above, once the relationship (in this case a linear relationship) has been estimated 
we can produce a linear equation in the following form: 



y — mx- 



(7.1) 



And once an analytic equation such as the one above has been determined, dependent variables at 
intermediate independent values can be computed. 



7.1.2 Performing Linear Regression 

Regression analysis with MATLAB is easy. The MATLAB Basic Fitting GUI allows us to inter- 
actively to do "curve fitting" which is a method to arrive at the best "straight line" fit for linear 
equations or the best curve fit for a polynomial up to the tenth degree. The procedure to perform a 
curve fitting with MATLAB is as follows: 

1 . Input the variables, 

2. Plot the data, 
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3. Initialize the Basic Fitting GUI, 

4. Select the desired regression analysis parameters. 

Example 7.1 

Using the data set above, determine the relationship between velocity and time. 

First, let us input the variables (Workspace > New variable) as shown in the following 
figures. 
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Figure 7.1: A new variable is created in the Workspace. 
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Figure 7.2: New variables are entered in the Variable Editor. 



Second, we will plot the data by typing in plot (time , velocity) at the MATLAB 
prompt. The following plot is generated, select Tools > Basic Fitting: 
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Figure 7.3: A plot is generated in Figure 1. The Basic Fitting tool can be initialized from Tools 
> Basic Fitting. 



In the "Basic Fitting" window, select "linear" and "Show equations". The best fitting linear 
line along with the corresponding equation are displayed on the plot: 
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Figure 7.4: Basic Fitting window is used to select the desired regression analysis parameters. 



Now let us do another curve fitting and obtain an equation for the function. Using that equation, 
we can evaluate the function at a desired value with polyval. 

Example 7.2 

The following is a collection of data for an iron-constantan thermocouple. 2 



2 Engineering Fundamentals and Problem Solving by Arvid R. Eide, Roland Jenison, Larry L. Northup, Steven K. 
Mikelson , McGraw-Hill Higher Education. ©2007 p.l 14 
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Temperature [C] 


Voltage [mV] 


50 


2.6 


100 


6.7 


150 


8.8 


200 


11.2 


300 


17.0 


400 


22.5 


500 


26 


600 


32.5 


700 


37.7 


800 


41 


900 


48 


1000 


55.2 



Table 7.2: Temperature [C] vs Voltage [mV] 

a. Plot a graph with Temperature as the independent variable. 

b. Determine the equation of the relationship using the Basic Fitting tools. 

c. Estimate the Voltage that corresponds to a Temperature of 650 C and 1 150 C. 

We will input the variables first: 
Temp= [50 ; 100 ; 150 ; 200 ; 300 ; 400 ; 500 ; 600 ; 700 ; 800 ; 900 ; 1000] 



Voltage=[2.6;6.7;8.8;11.2;17;22.5;36;32.5;37.7;41;48;55.2] 
To plot the graph, type in: 

plot (Temp, Voltage) 
We can now use the Plot Tools and Basic Fitting settings and determine the equation: 
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Figure 7.5: Basic Fitting window is used to select the desired regression analysis parameters. 



By clicking the right arrow twice at the bottom right corner on the Basic Fitting window, 
we can evaluate the function at a desired value. See the figure below which illustrates this 
process for the temperature value 1 150 C. 
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Figure 7.6: Estimating the Voltage that corresponds to a Temperature of 1 150 C. 



Now let us check our answer with a technique we learned earlier. As displayed on the plot, 
we have obtained the following equation: y = 0.053x+ 1.4 This equation can be entered 
as polynomial and evaluated at 650 and 1 150 as follows: 

> p= [0.053, 1.4] 



0.0530 1.4000 



> polyval(p,650) 



ans = 



35.8500 
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> polyval(p,1150) 



ans 



62.3500 



7.1.3 Summary of Key Points 

1 . Linear regression involves fitting the best straight line relationship to explain how the varia- 
tion in a dependent variable, y, depends on the variation in an independent variable, x, 

2. Basic Fitting GUI allows us to interactively perform curve fitting, 

3. Some of the plot fits available are linear, quadratic and cubic functions, 

4. Basic Fitting GUI can evaluate functions at given points. 



Chapter 8 

Publishing with MATLAB 



8.1 Generating Reports with MATLAB 1 




MATLAB includes an automatic report generator called publisher. The publisher publishes a script 
in several formats, including HTML, XML, MS Word and PowerPoint. The published file can 
contain the following: 



This content is available online at <http://cnx.org/content/m41457/Ll/>. 
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Commentary on the code, 

MATLAB code, 

Results of the executed code, including the Command Window output and figures created by 

the code. 



8.1.1 The publish Function 

The most basic syntax is publish('f ile' , 'format') where the m-file is called and executed 
line by line then saved to a file in specified format. All published files are placed in the html 
directory although the published output might be a doc file. 



8.1.2 Publishing with Editor 

The publisher is easily accessible from the Editor toolbar and file menu: 
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Figure 8.1: Publish button on the Editor toolbar 
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Figure 8.2: Publish item on the Editor file menu. 
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Example 8.1 

Write a simple script and publish it in an html file. 

Select File > New > Script to create an m-file. Once the editor is opened, type in the 
following code: 

x = [0:6]; % Create a row vector 
y = 1.6*x; °/„ Compute y as a function of x 
[x',y'] °/ Transpose vectors x and y 

plot (x,y), title ('Graph of y=f (x) ') ,xlabel( 'x') ,ylabel('f (x) ') ,grid °/„ Plot a graph 

Save the script as publishing. m and select File > Publish. An HTML file is generated as 
shown in the figure below: 
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Figure 8.3: A script published in html 



8.1.3 The Double Percentage °/„°/„ Sign 

The scripts sometimes can be very long and their readability might be reduced. To improve the 
publishing result, sections are introduced by adding descriptive lines to the script preceded by °/X 
Consider the following example. 
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Example 8.2 

Edit the script created in the example above to look like the code below: 

°/o°/o This file creates vectors, displays results and plots an x-y graph 
x = [0:6]; % Create a row vector 
y = 1.6*x; °/„ Compute y as a function of x 
•/.*/. Tabulated data 

[x',y'] °/ Transpose vectors x and y 

•/.•/. Graph of y=f(x) 
plot (x,y), title ('Graph of y=f (x) ') ,xlabel( 'x') ,ylabel('f (x) ') ,grid °/„ Plot a graph 

Save the script, a new HTML file is generated as shown in the figure below: 
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Figure 8.4: An html file with sections 



8.1.4 Summary of Key Points 

1. MATLAB can generate reports containing commentary on the code, MATLAB code and the 
results of the executed code, 
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2. The publisher generates a script in several formats, including HTML, XML, MS Word and 
PowerPoint. 

3. The Double Percentage °/ °/ can be used to creates hyper-linked sections. 
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Chapter 9 



Postscript 



Dear reader: 

Although I began working on this book awhile ago (c. September 2010), the current iteration of 
the manuscript is still far from being complete and I anticipate it will take a few years to bring 
it to a decent state. The text needs proof reading, extensive editing and more proof reading. I 
would therefore appreciate readers' suggestions, comments and sincere criticism emailed me at 
sbeyenir at my dot bcit dot ca. 
Thank you, 
Serhat Beyenir 



This content is available online at <http://cnx.org/content/m41658/Ll/>. 
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